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ABSTRACT 
The cavitation Derformance of a seven-Dbladed, highly 
skewed nrooelliar ee a Naval aux dot ak Y. (AQO-177>) is 
evaluatecac using Prec iri ng surface mumerical Hydteagymam ic 
methods. ATL el mpOGtant CCGMteritout ion re aqcurate predrit ton 


Se Gavitattion performance 1S Shown t¢90 Be an accurate model 


of the affective wake. 


A new fede: of the effective wake 1s described. A 
e000 am Leo calculate unsteady time-averaded, Owe Tot 
Cmeectimterentially averaged, face lod Bo int velocities 1s 
Dresented. The Droneller-induced velocities, along with 
tne Olie@ tt a tl nominal wake, are combined with Huang's 
axisymmetric effective wake scneme on ble-snaoned wake 


seaqments to determine the effective wake. 
The metnod is used to Stead ret severe C¢avitatton 


exgtent. This pradiction is confirmed by SSPA experiments. 
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i en ODUCT TON 


Ane David W. Tayviloct Naval Ship Research ang 


Development Center has recently designed a 
seven-Bladed, nLah y skewed propeller [Or a Naval 
Auxiliary CAO=]177). Tis design is instructive 
eoncerning the State of the act in computer-aided 


propeller desigqn. 


The Ni StOny 5Ot wene oropel ler G@esitanwis chronicled 
Piet Goof ts by Boswell (1), Valentine and Chase (23, and 
Hendrican and Remmers Car: Computer-aided propeller 
aopi tcGations Included Dieting Line theory, used to 
detarmine Piesmmecaagial Poadmdrs tf ibut 1 on and the radial 
hydrodynamic Deeen wang iva, ang rE ting sic tace theory. 
used to determine the final geometric Octet 


Sst rirbpution,; Gamber Gis trou tron, and final propeller 


offsets, Phe bua im Eto ive tes, trailing and leading edge 
details, 2040 1ona | thickness added to the Pre hi ing 
edge, tio geometry, rake, and nub details. Cavitat von 


performance was predicted Dy tne method of Baer tt io eane 


Emerson (4) and was also eavaluated by experiment. 


However, sea tttais of the AO-177 Pema awed a 
severe UNSteady Ppropetier torce proSlen. Inspection of 
the Dropelier indicated that 1k Was SGavitating 
signiziicantily worse than nag deen predicted. The 
Barret it and Emerson method, Surgqina lly designed oy) 





rour-to-six-SOSladed Dropelileérs with an expanded abed 
mae OO 8.60 ; aid fre. extend well to thea saven-bladed 


AO-177 pronelier with an expanded abies) fa ‘to voy. 8 ey cere 


even though both were POG merchant-type hulis. 
Furtner, scale eaeaftacts could have thrown ort the model 
tests. 





je a Pe A Gowme ress EEFECTIVE WAKE 


The AO-177 propeller performance was avaluated at 


ge ae & using some humectcical propeller analysis computer 


programs developed at MIT. Kerwin and Lee (5) describe 
MiT-PUF-2, a computer program to Galcutate forces, 
steady ang unsteady, generated by a propeller. 


Miter rv; a steady field point velocity po Og ram, i's 


GescErhped 1m a2 report by Iya ps g C6). Leas C7) freports on 
Mer—PUE=3, 2 computer program to Predict steady and 
unsteady Bee oO piece l et Caviwatl Onewsextemt “and volume. All 
three boOGgrams « mode i the propeller as = aqrid Of 
discrete vortex seqments. The wake is also modelled as 


a vortex grid tin each program. 


Gavitation performance Ga leoubet 1ons , more tn am 
Propeller FOLGE and tnduced velocity Galeulat tons, are 
highly sensitive to the wake data used as tn pce 
Experience at 1g OO a and elsewhere has shown that using 
the nominal wake data to approximate Phil ow COnG1 tions , 
while acceptable for Sf ooel ler rorce and induced 
memoecity™ caleulations, is mot acceptable for cavitation 
performance calculations. The wake velocities measured 
Behind a model in the absence Of any Dropeller (C(thea 


nmominal wake) Must | demoed ted CoO account for change in 


Boundary Eonartions imposed by the presence Ot a 
Propeller. This new, modified wake is generaliy called 
the “etfective wake. ' There is, unfortunately, no 





general method EO: atic. © lpa. Lilg the etEtective wake 
generated By a Dar Creu lbar Pprooeller mi a general 


mominal wake. 


Hitandusanad OtChers s<a), 7) have developed porn theosy 
and mlumerical schemes tor soa leulat ing tre effective 
wake of a propeller im an axisymmetric nominal wake. 
Deteon cote Veritas fas attempted to adaot Huang's 
method to non-axkisvmmetric wakes, such as surface shigpn 
wakes, Biv Gi val din ¢ the nominal wake into pie-snhaoped 
segments and then applying the axisymmetric wake 
Salculationmwithin eacn pie se@aqament, using tne nominal 
wake and the steady propeller induced velocities. MaAG 
aoproach is an improvement over blindly apolying the 
akisymmetric etfective wake Galcul aiwion to a surface 
ship's waxes, But, owing to the very sharp AO-177 wake 
(aven DY Surface ship standards), it was judged that 
some further refinement Of Huang's metnaod was 


Fequtired. 


Huang's method assumed an axisymmetric nominal 
wake and a Sh each propellar-induced velocity tiela. 
The Det Norske Veritas method assumed a general nominal 
wake and a steady propeller-induced velocity Eel a 
Ein ct of a completely general effective wake 
Calculation, the most obvious refinement was to assume 
a general nominal wake, but allow an unsteady 


Pproveller-induced Velocity t1eld. This was the method 





Woeeiag LOL tEhis paper. 


he current method of effective wake caleulations 
can be summarized in the following steps: 

Lo Use the nominal wake ang oroopelier data to 
calculate the unsteady propeller-induced velocity field 
at representative 0 Ojimpcoee to oe fi inside and outside the 
propeller disk. 

Zo Divide the mominal wake into small pie-shaped 
Sseqments as in the Det Norske Veritas method. 

<)), Calculate the effective wake velocities within 
each pie segment using the nominal wake velocities and 
the unsteady propeller-induced velocities, assuming 
that each pie seaqment acts as though it were part of an 
akisymmetric wake with the same Ea dii-at WellO G6 It y 


Gpstribution. 





Ge ea MElHOD OF EFFECTIVE WAKE CALCULATION 


The two new aspects Of the currant effactive wake 
@eiculations are tne propealler-induced unsteady fireld point 
velocity CUE PV. Gaul CU la telom and the use of these 
time-averaged But not circumftferentiallv-averaged velocities 
in a pie segment nominal wake modification Ga toutat 1.67 
So LEWAKE). Eacn of these two new methods wili bea 


Mimmioct tated by a typical calculation. 


Tne unsteady ELelea ,oorne velocity CUP RV > (.p f od tam Ls 


bastcally a Generaitzation of the Steady field paint 


welocity Drogram (MIT-FPV) bY Ib rogl (0 oe The unsteady 
version uses the fu lel wake lattice arrangement also 
emploved in MIT~PUF-2. PErOdtam 1hDUtS are the FILE14 DATA 


file output from MIT-PUF-2 and interactive Girections : 
Program outputs are either Single Dorn t Calculations 
directed to the user's terminal or ardata £11 to be used 


momeinput into PIEWAKE, for arfective wake calculations. 


UFPV first “reads the FILE1¢ DATA input, containing 
mostly Propelier geometry, Singular: ty strengths, and 
Singularity geometry. Ptw should be mooted that the number 


et propeller revolution tima steps specified in MItT+-PUF-2 
must be divisible by the number of propeller blades . For 
eftective wake Galceculations, the wake tnbut to MIT=-PUE-2 
Srgant Strictly to Be the effective wake; however, using the 


mominal wake to calculata the pPpropelier-induced velocities 





mama GOOQd abDOrogkgimation. 


The user next a elia, (ele <j Ered gd Hemme COOnTdinates. The 
field point position may be entered Big evaiitmar teat 
coordinates or in terms of theo fotating oropeller and wake 


ar ids. 

The velocity induced By one Blade at Baer time step is 
then Galéulated. At each time steo, the vortex ITattice is 
assigned the strengths corresponding to the time step of 
the propeller Ee wat ton. The velocity imduced DY eaca 
eiement of the vortex GG iG, born IO lade: and wake, is then 
summed to obtain the total velocity induced by one Db iaae ini 


<n a © particular angular position. 


Mmemvelocitires induced by each Blade of the propelier 
are tne) $uTamed to produce the torts propeller-induced 
Verlocity for that Bosition of "tne propeller. Velocities 
are Calculated for diffarent propeller positions and then 
macmonicaliy analyzed. The Zenmoetin Raprmonic, them, is the 
time-averaged Bropetletr—-induced velocity for that field 
point. (By eer all points at a given radius we could 
obtain the Circumferentirally-averaged, time-averaged, 
propeller-induced velocity, the Oto ut Qt the steady 


uy a a 


For effective wake calculations, such velocities are 
computed for up to -60 points par radius for Radi Varying 


Prom propeller hub to L Soy times the propeller radius. 





Between the PLOpel ter Nupeand propeller tip, the velocities 
are eatrewlated at the leadina edge panel on the pnropeiler 


geld. 


For the AQ-177 propeller analysis, 36 points per 
radius for eiaht radit were used to St Gat Gara, £ 1te of 


unsteady field point velocittes. 


PepwAKRE iS 4 program to perform Huang's axrisymmetric 
effective wake contraction calculations hg eacn of many 
Ple-snhaped segments of the wake. rat uses nominal wake aata 
and unsteady propeller-induced Veloe1<.y gata trom Urey. as 
MPaeuts to sroduce an @rrective wake valocrrty data fiie as 
output. 

After reading nominal wake data and the induced 
weroctty Prelag data, PILEWAME Teztrapolates the mominal wake 
data to the Rub. The imnermost radius of the induced 
velocity field is taken as the Aub. (UFPYV had taken, as 
ites inmeermost radius, the RUB "Tacgius used in MiIT=-PUF-2z and 
passed to UFPV Dy way Of Thee lo -dgats f£ fie. 2 

oA 

in each oe seqment, Huang's method is used to 
calculate an etfectiva velocity and themetrective radius 
corresponding to that effective velocity. if “tre Nwinéeésr ica i 


orem, this method is tterative in nature. 


A Trin i te differance equation, given in reference (9), 


is used to calculate the effective welocnt ty . Given tne 





mominal wake at radit Gacstimed= of sthe fEirst tterattion, to 
be the nominal radii) in the presence of aepropes |r, an 
ettrective Vekoc lt ¥ Gan De calculated Et an errective 
velocity at the next radius away from the hub is known. 
The eftective velocity at wit hmemrowt ormos t radius is assumed 
ae) be equal to the mominal velocity there, enabling all 


etfective velocities to be calculated. 


A second finite difference equation, also qiven ig! 
Treterence C2), is used to calculate the effective radius 


corresponding to the set of effective velocities just 


eeoerettlated. Timi s equation uses nmominal Badivi, nominal! 
velocities, and apparent velocities. (An apparent velocity 
is the sum Of the- effective vee ore 1tly and the 
propeller-induced veioci ty at a DO Lit 22) Given one 


errective radius, the Meet Eradius away from the nub can be 
Seevctulated. Sinc¢ee the iMmermost cadius On aa-l velocity 
tields is set au vGMe hNuo wD aaiwts,nraklt @ffective radit can be 
determined Secuettiatly from inmermost to outermost. These 
Mmeori Can tnen be used in the finite difference equation to 
Sea iculate the mext Lterat ion of effective velocities. 


PIEWAKE iterates fiva times for each pie segment. 


in addi1ti0onm to thea effective velocity wake flerid:, 
PIEWAKE also computes the mean YVeboctrty, tne volumetric 
mean velocity of the mominal wake, the volumetric mean 


velocity ofr the effective wake, and the effective blockage. 


The effective velocity field then then be processed 





ema  usedgd by MIT=-PUF-3 as any other waka field. 





IV. AO eet oom raewmm se CALCULATIONS 


Appendix l Contains plots [= abe mMominal, induced, and 
eftrective velocities varying Elrcwumicremt ial lyv-cat Saveral 
specitied Prac Appendix jk Contains the same 
Pmiotrmation, exceot that the nominal and etractive 


welrocrtres aré normalized on the respective volumetric mean 


velocity. These plots COnMtTirm that Ehe moditications £10 
the nominal wake outlined gg Chapter Pict are indeed 
reasonable. Several saPrent conclusions can be drawn from 


these plots: 
1) Er hect twee veloe it y Modriications of the nominal 
velocities are in the expected airect ion. Physically, tna 


Pmeosition of the propeller on the nominal wake draws more 


water into the slipstream, thus increasing tne velocities. 
TAs "faster water" will come from PagGgii “Greater than the 
bield point being sampled; im -.otnern Words, the wake will 


Semem to “contract.” 


Z) Effective veloc ty ModLtrtercat ons of the nominal 
wake velocities are Greatest at lower radil. The wake 
a 


meomtraction”™” will seem to concentrate here. 

3) In general, the waxe peax is more narrow in tna 
effective wake than in Cie nomirna ft wake. in. “ene Gace of 
the A0-177 nominal wake, the existence in the presence oft 4a 
propeller of a waka peak as sharp as that in the nominal 
wake is counter-intuitive. Tier etore:, a decrease hogl the 
wake peak is to be expected. 


4> The normalized effective Velocities at 0 degrees 





are moto ep Dreciably difterent fc om that of the mMominat 


wake. This is disappotnting and was mot erpected . [t may 
be a result of © arceieaut wor tt Cult y ; LTenored DY tone 
axisymmetric erfective wake CayvecuUlawuron pmertotmed in each 
pie segment. This assumption is least approoriate in the 0 
degree region or the wake. A more general efteaective wake 


calculation scheme may be fequired to obtain better results 


Pie this region. 


Results Ort Mit-PUF=3s “cavitation performance of the 
AO-177 propeller are shown in Apoendices [IT and IV. AS an 
Prrrcator of the influence of effective wake modifications 
to the nominal wake, MI{T-PUF-3 was run twice, once with the 
mominal wake as ingout and once with the effective wake (the 
more ageieurate boumdary condition > [2e input. Appendix Til 
Shows cavitation performance elo Gi—g! the nominal wake as 
Poe 12 t while Apoendix IV snows Savit ation performance “diven 
the effective wak@ as input. Imescgenertcal, the plots show 


taiat tire cavitation extent is not changed much, aiLthouda 


cavitation performance is better tin the effective wake at 


the lower radii wnere, as has Been shown, the erfeactive 
wake moms Pee howe are the Greatest. AtomMme ,vouter Gadus 
where the effective wake Mogdytreat ion is the least 
peeamiticant, the Carey it tacron Deb ronrmance is least changed 


between the nominal wake and effective wake cases. 


Appendig V contains a er aon comparing Gavitation 
Sxtemt of MIT-PUF-3 output given nominal wake tinbut against 


that given effective wake input. Tt plots the length of 





SGavitation on the Seventh radial Saviteation oOanel Cabout 
the .8 propeller radius} mon-dimensionalized am the loc at 
Sa1erd length as the blade rotates fmt ougqnh 340 degrees. 
ais aqraon Indicates that the effective wake modification 
did Meceneduce cavitation extent on a Dlade, a result noted 
SeELier. However, 4a larger contrast Between the nominal 
wake and effective wake cases is shown in Appendiz VI which 
@eagnhs the cavitatron volume oni a blade in the two cases. 
The volume is non-dimensionalized on the cube Gre the 
propeller radius. TA i.s grapn shows a droo of apnrozimately 
35% o £ peak Gavi at ton volume. hrs change rasa Gavity 
volume Gemonstrates the sensitivity of Cav tte t iON 
performance to small changes in the wake profile. ThAis 
Somcirms that ettorts to refine etfective wake calculations 


were indeed justified. 


The effective wake cavitation performance indicated by 
Mi{T-PUF-3 compares tTeasonably well with experiments 
performed By tthe SSPA propeller tunnel. These experiments 


Measured cavitation performance of the AQO-177 prooelier Dy 


Photographing model tests. While no cavitation volume 
2 

results were obtained, CavrpeatLronm extent. could De compared 

with M1IT-PUF-3 pradictions. Despite reductions Of the 


effective wake S2aiculations compared with nmominal wake 
momreulations, MIT-PUF-3 still tended to Htrediect Greater 


Cavitation than the experiments indicated, particularly at 


the middle Bagi: DaLs seems to in dt6 4.2.2 that the 
etrfective wake modifications, while tending to the proper 
direction, were not large anough. This tS consistent with 





the Gonelusions from comparing the two wakes that a2 more 
general effective wake Schemaemerimec  iding cCadial vorticity in 


high wake regions is required. 


\ 





ae CONCLUSIONS 


1) Gauge 7 iat io nh of an ettective wake using Huana's 
aKkRisymmetric method el pie-shaped segment with unsteady 
propeller-induced velocities is a reasonable approximation. 


rt er ill tands to underestimate the effactive wake 


Mmoaitication Be the nominai wake Lon the Vie 1 Mi-c y of 9 
degrees. eomcOrrec t tHilswe aemore general effective wake 


eartculation Scnemew ine icing faqatral vortterty etrects may 


be redauired. 


9 The current effective wake Ca leulation scheme can 
De used to predict cavitation performance tor Aighiy skewed 


mropellers with up to seven blades. 


3? Cavitation performance is indeed sensitive to 
small changes FS pl the wake Or of ite. Seemingly small 
improvements in eftective wake Galeulations yield 
Beomifticantly more accurate Gavi ta6 1.0n performance 


Beedictions. 





Vols RECOMMENDATIONS 


1) A more general effective wake calculation scheme, 


Meremrerdingd the errects of radial vortictty, its required. 


Ze The computational time Camad “rence "cos t ? of the 
unsteady tield Botnet Velocrtey program. (UFFV > can orobably 
be cut without Sram if icane sacrifice Or accutac y By 


decreasing the RuUumMper of points Satleulated per Paqgiwus . 


Appendix I shows that the vVarbisacion Of induced velocity 
Been rotation ts Solita toad smooth. Ea Cew laste bon., Of fewer 
soints and estimattion of the Femalnindg points using 


MasmMOnic analysts would probably not reduce accuracy. 





RX APPEND EX fo 5% 


Nomimal. Induced, and Effective Velocities 





° 


ar ©. | 







TOL EFELE RTL 
Sonn (heen 
CONE 
COCA Ve 
COO CSSERZCCe 
CEES Per 
SSSRREDD2 92 
COO 
COO ae 


BREED AREER 
\ SERRRE IEEE 












CTOAXAS CEDIA SAYA WWIXY 


v( VOL) 


Nom 1.0 
1,0 
{.0 


Ind 
Eff 


= 


.4@ RADIUS 


NOK LND. 7 EPP VELOCITIES, 













Lit terry |} NET 
CECH Eer REE 
SRR 
Hat Ye 






Sap 2aknn 
COO CCC 
COCA 
CCCCSSSECr aT 


WN 
Q Q & & Q a>] © & Le ] © 
I t 









CTOADSASCYIA = SAYM TWIXY 


v(VOL) 


Nom 1,0 
Ind 1,0 
Ber 220 


~ 23 - 


.68 RADIUS 


NOM, IND, EFF VELOCITIES, 





90 





JERR EREREERIGEE 
LPL i eet et 





CIGCAXA/SCYDA =SAVN IVIXY 


wm 

r= 

© 

= 2 2 © 

~~! € 8 8 

— wt ed 
Ew w 
OS & 
= iH 


a Oi 


.78 RADIUS 


IND, EFF VELOCITIES, 


NOM, 








COC 
CECA 
COCR ETT 


” N _ a o) roe) ms oO 
_ _- ~ = A) ‘@ 2 ~) ) A) Q rm) re) A) e 





CTOADSASCYIA =SXAYM TWIXY 


~o_-_ 

— 

© 

= Oo = © 

= 68 8 «@ 

- SS 
Ed Ss 
O & “ 
= ei (x) 


= 2g = 


-Q, 


-150 -120 -90 


~189 


.80 RADIUS 


err VELOCITIES: 


IND, 


NOM, 





160 














Tlie te AT 
FEEAEEEEE EEC 
Tet EE EE ET 
See eee 


68 80 {20 1S 


38 


~30 


-60 





ria 
BEER 
CCC 
CCE 
LOOSE 


z § ] 


-~{SQO -120 -98 


480 


CTOADAS CYDIA = SAYA IWIXY 


v(VOL) 


Nom 1,0 
1.0 


Ind 
Eff 1,0 


’ 
tO 
ON 

a 


.98 RADIUS. 


EFF VELOCITIES, 


IND, 


NOM, 





Tee 
S00) eee 
CCA 
CCEA 
COCO SSECCC Ee 
COCR 
COO EE 
COCA Ce 
CECA CCC 
CCA 
SCOPE 
COO EE 


= = — — & fe Q ® & a] Q & & & 


180 





152 





120 






92 






6 






38 












-30 


-68 






-98 






-120 


-150 


= 


-{80 


-0.2 


CTOADA/CHIA = BAVYM TWWIXY 


OO 


a e 
rt =e 


v(VOL) 
1,0 


TH = 


Ma 


27 


Nom 
n 


.98 RADIUS 


NOM, IND, EFF VELOCITIES, 





Nominal, 


AAS PEND Xe ox x 


PmMcwced, a4. 2: tect iva 


Non-Dimensionalizad 


\ 


Velocities, 








Epe aes Nee 
a ( coc 
Stitt TT ET Tt 
BENSSEREEE Ae 
PEL PS ARNZE ETT 
Pitt dy | Peal 
SERRE EDEE Scene 
LET rer IN | Td 
Hey | a - RE 


a) PT UT TAT Ty 
ENE ERR 
ime 


_ = = —_ o 


152 









120 














-90 






~1{80 -15@ ~-122 


=-@.2 


CIOADA/CADIA = 3NAVM IVIXY 


V(VOL) 
Nom 0.732 


Ind 1,0 
ber 0,787 


Sis) 


180 


.48 RADIUS 


NOM, IND, EFF VELOCITIES, 








RCCCCCC 
USSSSSRSRR 
PPRBEE TT 


-~ _ — ~ far) © © o o © @ © © ° o © 
i 





© 
© 


CIOASASCUIA §=SAYM TWIXY 


v(VOL) 
Nom 0.732 
1,0 


Ind 
eff 0,787 


- 30 - 


,6@ RAOIYS 


EFF VELOCITIES, 


IND, 


NOM, 









ARE CEEEEE EE ee 
RECEP EEE 
NGS ve 
CPSSEECT EY PEL 
COCA SSL 
CCE eer TC 
CBSE SEL 
CACC ECL 
STEELE 


TET tT 
We 
Prostate | | | Le] dd 


sa = = == ° @ @ o o @ @ @ Q o ° far 








690 







38 











CTOAXA/S CYDIA = SAYM TYIXY 


V(VOL) 
Nom 0.732 


Ind 1.0 
Eff 0,787 


i 
Lor] 
raair¢ 

§ 


.72 RADIUS 


NOM, IND, EFF VELOCITIES, 









AME EE PPE ey 
Seb} 





AAPL TT ET EEE 
OPT TT TTT ET tt 
fe 


- - - - 9 Ea rar) Q rar) 9 © @ rer) rer) S 






-Q 


CTIOADASCSNDIA = SAYM IVIXY 


v(VOL) 


Nom 0.732 
1.0 
0.787 


Ind 
Eff 


oo 


LG ark 


99 120 1659 189 


68 


39 


-120 -90 -60 #£-39 


-150 


.829 RADIUS 


NOM, INO, EFF VELOCITIES, 





98 120 1S9 168 


69 






38 






A 


CCE Ere 
aes 
RICCO 
A 
pop 


= 30 


~69 









-99 







> 


- - — - © far) Q >] Y oS fe ] fer] pa] a*] - far) 







CTOADASCYIA BAYA TWIXKY 


F oY 

mi A ~ 
Ov” ee) 
-~ ~~ O 0. 

~=—_—_Z" © ¢ @ 
—- OwnA Od 
Ema ws 
Of WwW 
=e & 


4 
LS) 
Ww 

§ 


EFF VELOCITIES, .9@ RADIUS 


NOM, IND, 





188 


G2RRRRRRRE 
DAE EE EE Tt yy 






98 


30 







CPT ECC 
Cet tT tT EE ET EE 
dV et tt ETT TT TT Ty 
PPE PET Ett ttt yy 
> SR 


-~69 







-92 





-122 





-{S2 





Fe 


! 
1.8 
8.9 
Q 
Q.7 
Q 
Q 
Q.4 
Q 
@.2 
Q 
2.8 
Q 
-2 


CTOASAS CYDIA = =SAYM TVYIXY 


V(VOL) 
Nom 0.732 


Ind 1.0 
,Eff 0.787 


= 


15@ 


122 


69 


.98 RADIUS 


-38 


NOM, IND, EFF VELOCITIES, 





MA TAPPENDIA L122 ** 


Piro e 116 t Cavitation, Nominal Wake 























A ee 


MT \XIN A\ 


= 
\ I Wa<s 
or) 


( 
ASL 


\\ WX OF 
WG 








AN | 
RIN 5 
oo 

x <A 





ccs 
SS 


| 
\ 
















Zi 


iit 
aN | a \ 


XS A 








aR 








eae rer ENDL x Ly ex 


PrOug@ulLerebavitation, -Etfective Wake 


Say 








Tay 

\\ £- 
PN a 
S7] 
_ —" % 


. 
. 
» 

















MO & 
ee 


\ 
/, 
@ 
N 
4 CW 
4 r= 
A © 
we! 
Nl Zs — 
w=) 















Rik 











Ue \% ON 
NO \ \s 
KO 











eA APE OND A Vi ex 


Savirvebemagtnm, Nominal and Eftective Wakes 








ep Win) av\b LOS isle)" Sys V/NGOSIO)N 
BUSINES als 


OL¢ Obe CRRA 08 | Oa CGO Oe) 0 OC-10e] Oe. 


+ FREM SS Y) 


A phi 


bad we ee - 


Mevlunte ClbsINaoy finer seeaiae 


. | Sa 

| ace | : 
PS BAP POAT TA. [a 8 

oe f r = 


us . pe dinaeae): BQ" 
92H] - ; 
teury ON ‘ | 
| = 
9: 





<x SAE ENGR yl x 


GavrtyeVoLuUme, Nominal and Effective Wakes 








SPW ie eel] a Se VN LION 
BSN Say 1S 


BLZ Ord BIZ O81 BSI OZ} O6 OS BF BOB OE- BO- AG- 


_ * PASSE BSCE pes ebeneeeepesesene ass” \——- ABOVE Peseses Z1Z1cm 





JUEGO SPN SIGN 7S IIe es) fike 


a 200 

a bOO’ 

— - 

See es Ist 200° 

* ar B10" 
= cal , 210 

|| ——-}-— bi 

é — _ ne ons 

_ —- 810° 


— EXHAU/SA 


- 56 = 





[noe END AV lok % 


Sorerneesc perimental FPropelier Cavitation Extent Piots 

















eX APPENDIA Viil *x 


Nominal Wake Diagram 








tt deere tT 
SH 












COOH SSE Cr 
CCCP 
CCC EBEET 
an) 49 7 SR00e0ne 
CCA or 
EE E 
















CTOADAS CYDIA SAYA IWIXY 


- 61 «- 


= CC 
CCS 


V(VOL)=1.0) 


Actual Velocities (I.e., 


NOMINAL WAKE, .4, .6, .8, .98 RADII 





xx APPENDIX 14 


Etfective Wake Diagram 



















li i@at itt 
CCH «*-2 50000 
iN pti ttt ty 
SERA 
tt] RSS 
Pid dT Ss. 


- “SEER 
eae 





® 

Pa 3) 

te | | | 

ia © 

fe) 

i 

&) 

Oo 

} 

fa] 

WN 

| 

® 

Fe) 

| 

& 

ie @) 

2] ro | _ ®& oO) bee] Na 
° i 


CTOADA/CYIA = SAVM IVIXY 


- 63- 


Actual Velocities (I.e., V(VOL)=1.0) 


8, .98 RADII 


63 


4, 


EFFECTIVE WAKE, 





xe APPENDIX XX *% 


UES vVeErogram b&1S5t ing 








NQNAIANAAANANAANAAa 


aan an00 an0 aan aan 


aAaaqnNn 


aan 


“4 


- 65- 


ee Reanakhenaeanennweaneee eee eee ees ee enw staneewnaaeeenanaeeee UF POOO10 

* = UF POOO20 

* MAIN PROGRAM IN FPUV, A PROGRAM TO CALCULATE * UF PQOO30 

* UNSTEAOY VELOCITY HARMONICS AT FIELO POINTS bg UF POOO040 

= SPECIFIEO INTERACTIVELY IN ROTERM. * UFPOOOSO 

* * UF POOO6O 

ee ee a ee ee a ee ee ee a ae UF POOCO7O 

UF POOO80O 

UFPOQOOSO 

**« COMMON TO CONTAIN GEOMETRIC PARAMETERS UF POO 100 
UFPOO110 

COMMON / GEOCOM / MM,MN,NN,NW,NTORL,X(12,10),.Y(12,10),2(12,10), UFPOO120 

& AUYeeton 1 ij y rire ( 10, 11), ZTIF( 10,511). XH( 101.7). YHC 101, 7), UF POO130 

& Zit 401.7), XW02 1), 10) , YWC 21, 10) .zw(21.10).RH UFPOO140 

UF POO 150 

*#* COMMON TO CONTAIN SINGULARITY STRENGTHS UF POO 160 
UFPOO170 

COMMON / SINCOM / GT(10,10,60),68(10,9,60),GTW(10),S8(10,9),GTV UF POO 180 

& .,GAM(9,60),GMEAN(9) , OKFW(20) ,HBWX ,HBWY ,HBWZ UF POO 1390 

UF POO200 

*##*« COMMON TO CONTAIN FIELO POINT OATA, NBLAOE, IOENT, NSR, ANO VSR UFPOO210 
UF POO220 

COMMON / FPCOM / N&LAOE,NSR,XFP,RFEP,TZFP, YBAR,ZBAR, UF PQO230 

& OBLAOE ,OTFP ,»IOGENT(18),.TFP(60),VSR,IPLOT,OIAM UFPOO240 

UF POO250 

*** COMMON TO CONTAIN VELOCITY OATA UF POQO260 
° UF POO270 

COMMON / VELCOM / VIX(60),VIT(60),VIR(60),U(60,3),.NSAMP UFPOO280 

UF POQO290 

zee COMMON TO CONTAIN HARMONICS COEFFICIENTS UF POQO300 
UF POOS 10 

COMMON / HARCOM / A(15,3),8(15,3),AMP(15,.3),PH(15,3).NH UF POO320 
UFPOO330 

**#m TFTSTOP IS A FLAG SET INTERACTIVELY IN ROTERM TO INOICATE ENO OF PROUFPOO340 
; UF POO3SO0 

ISTOP=0 UF POO360 
UFPCO370 

*e* ROFILE WILL REAO OISK OATA FILE CREATEO BY PUF-2. NOTE THAT NSR IUFPOO380 
PUF-2 INPUT MUST BE DIVISIBLE BY THE NUMBER OF BLAOES, NBLAOE. UF POQO3390 

UF PQOO400 

CALL ROFILE UF POO4 10 

UF POO420 

wee INTERACTIVELY REAO PARAMETERS FOR TERMINAL SESSION UF POO430 
UF POO440 

90S FORMAT(/’ ENTER NBLAOE, NTORL, IPLOT.’/ UF POO450 
& ’ FOR HELP, ENTER ZERO FOR NBLAOCE.’” ) UFPOO460 

30 WRITE(6,905) UFPOQO470 
REAOQ(S,*) NBLAOE,NTORL,IPLOT UFPOO480 

IF (NBLAGE.LT.O) GOTO 20 UF PO04390 

IF (NBLAOE.EQ.0) GOTO 30 UF POOSOO 

IF (IPLOT.EQ.2) CALL PIEOAT(ISTOP) UFPOOS10 

IF (ISTOP.EQ.1) GOTO 20 UE FCOs 20 

UF POOS30 
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a00 
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aqanan 


s** ROTERM GETS OATA INTERACTIVELY FROM TERMINAL ANO COMPUTES UFPOO540 
NECESSARY FIELD POINT GEOMETRY. UFPOOS50 

UF POOS6O 

10 CALL ROTERM(ISTOP) UFPOOS570 
UFPOOS80 

*e* IF OPERATOR INDICATES ENO OF PROGRAM, ISTOP WILL BE SET TO 1 UFPOOS5390 
UFPOO600 

IF (I1STOP.GT.O) GOTO 20 UF POO6 10 
UFPOO620 

*** RUN THROUGH FPSTEP FOR EACH “TIME STEP.” UF POO630 
UFPOO640 

0D 100 K=1,NSR UF POO650 

UF POO660 

**2 FPSTEP CALCULATED VELOCITY INOUCED AT FIELO POINT BY ONE BLADE UF POO670 
AT “TIME STEP®* K UF POO680 
UFPOO630 

100 CALL FPSTEP(K) UFPOO7OO 
UFPOO710 

axe SUMVEL WILL SUM FIELO POINT VELOCITIES INOUCED BY EACH BLADE. UFPOO720 
UF POO730 

CALL SUMVEL UFPOO740 
UFPOO750 

**= FOUCAL WILL GENERATE HARMONIC SOMES NTs FOR INOUCED FIELO POINTUFPOO76O 
VELOCITY AS PROPELLER ROTATES. UFPOO770 
UFPOO780 

CALL FOUCAL UFPOO730 
UFPOO800 

*** PLOTVL PLOTS INOUCED FIELD POINT VELOCITIES AS A FUNCTION OF BLADEUFPOOS8 10 
ROTATION. UFPOO820 

UF POO830 

SO1 FORMAT(’ AXIAL CODROINATE’ ,4X,F5.3/% RADIAL COORDINATE’ ,2X,F6.3/ UFPOO840 
&® %’ ANGULAR COGROINATE’ ,2X,F5.1///’ AXIAL HARMONIC’ ,6X,F8.5/ UFPOO850 

& % RADIAL HARMONIC’ ,5X,F8.5/’% TANGENTIAL HARMONIC ’,F8.5) UFPOO860 
WRITE(6,901) XFP,RFP,TZFP, (AMP(1,JU),.J=1,3) UF POO870 

IF (IPLOT.EQ.1) CALL PLOTVL(NSLADE,NSR,NH,TFP,TZFP,U,A,8,AMP,PH, UFPOO880 

& XFP,RFEP, IDENT) UFPOO830 
GOTO 10 UFPOO3S00 

906 FORMAT(/’ NBLADE < O => STOP’ UFPOOS910 
& FOX: ‘=z O => HELP’/8X,’> O => NUMBER OF BLADES’/ UFPOOS20 

& ’ NTORL = O => TOTAL EFFECT’/8X,’= 1 => THICKNESS ONLY’/ UF POOS30 

& 8X,’= 2 => LOADING ONLY’/’ IPLOT = O => DATA TO TERMINAL’? / UFPOOS4O 

& 8X,’2 1 => PLOT -- NOT IMPLEMENTEO 05/02/82 -- RLU’ UF POOS5O 

& /8X,’= 2 => CREATE PIEWAKE DATA FILE’/) UFPOOS6O 

SO WRITE(6,906) UFPOO970 
GOTO 30 UFPOOS80 

20 sToP UFPOOSSO 
END UFPO1000 
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SUBROUTINE FPVEL(K AF .YF.ZF.VIX,VIY.VIZ) 
**« COMMON TO CONTAIN GEOMETRIC PARAMETERS 


COMMON / GECCOM / MM,MN,NN,NW,NTORL,X(12,10),Y(12,10),2(12,10), 
& <preeiO, 11).Y1T1P( 10, 11), ZTIPC 10. 11) .XHO101.7), YHO101,7)., 
& ZHiiO4 7). XW(21, 10). YW( 21. 10).Z2W(21,10).RH 


**« COMMON TO CONTAIN SINGULARITY STRENGTHS 


COMMON / SINCOM / GT(10,10,60),G68(10,9,60),GTW(10),$8(10,9),GTV 
& ,GAM(9,60),GMEAN(9),OKFW(20) ,HBWX,HEWY ,HEWZ , 


wee COMMON TO CONTAIN FIELO POINT CATA, NBLAOE, IOQENT, NSR, ANO VSR 


COMMON / FPCOM / NBLAOE,NSR,XFP,RFP,TZFP,YBAR, ZBAR, 
&  OBLAOE,OTFP,IOENT(18),TFP(60),VSR,IPLOT,OIAM 


DIMENSION Utx(20, 2) ,UTY( 20, 2) .UTZ(20, 2) ,UBx( 20) ,UBY( 20) ,UuBZ( 20) 
& |, VOFWX(20,9),VOFWY(20,9),VOFWZ(20,9) 
wee INITIALIZE VARIABLES 
MIO=MN/2 
NVFW=NW 
UMPFW= 1 
WAKX=0. 
WAKY=0O. 
WAKZ=O. 
VIX=0.0 
VIY=0.0 
VIZ=0.0 
DO 1 M=1,MM 
CO 1 N=1,NN 


----- CHOROWISE VORTICES ON BLAQE-------------------------------------- 


IF(NTORL.EQ.1) GO TO 8 
CALL VORSEG(XF,YF,ZF,X(N,M),Y(N,M),Z(N,M),XO(N+1,M),YCN+1,M), 
1 2ZOCN4#1,M),CVX,CVY,CVZ,SX,SY,SZ,0) 


----SPANWISE VORTICES ANO SOURCES ON BLAQE--------------------------- 


CALL VORSEG(XF ,YF.ZF,X(N,M),YC(N,M),ZON,M).X(N,Mt1),YON,M41), 
1 Z(N,M+1),SVX,SVY,SVZ,SX,SY,S$Z,1) 
SBNM=SB(N,M) 
GBNM=GB(N,M,K) 
GTNM=GT(N,M,K) 
IF(NTORL.EQO.0O) GO TO 9 
IF(NTORL.EQ.1) GBNM=0.0 
IF(NTORL.EQ.1) GTNM=0.0 
IF(NTORL.EQ.2) SBNM=0.0 
VIX2VIX+GBNM*SVX+SBNM*SX+GTNM*CVX 
VIY2VIY+GBNM*SVY+SBNM*SY+GTNM*CVY 
VIZ=*VIZ+GBNM*SVZ+SBNM#*SZ+GTNM#*CVZ 
CONTINUE 
----SEPARATEO TIP VORTEX ON BLAQE------ 2-7 ene n ne ene ree reer ee 
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UFPO1010 
UFPO1020 
UFPO1030 
UFPO1040 
UF PO1050 
UFPOQ1060 
UFPOQ1070 
UF PO1080 
UFPQ1090 
UFPO1 100 
UFO 1 16 
UFPO1 120 
UFPO1130 
UFPO1140 
UFPO1 150 
UFPO1160 
UFPO1170 
UFPO1180 
UFPO1190 
UFPOQ1200 
UErPO12Z10 
UFPO1220 
UFPO1230 


“UFPO1240 


UpPOtZs0 
UEPO 260 
UFPO1270 
UFPO1280 
UFPOQ1290 
UFPO1300 
UFPO1310 
UFPO1320 
UFPO1330 
UFPO1340 
UFPO1350 
UFPO1360 
UFFPO 1370 
UFPO1380 
UFPO1390 
UFPO1400 
UFPO1410 
UFPO1420 
UFPQO1430 
UFPO1440 
UFPO1450 
UF PO1460 
UFPO1470 
UFPO1480 
UFPO1490 
UFPO1S00 
UFPO1510 
UF PO1S20 
UFPOQO1530 





Cc 


IF(NTORL.£O.1) GO TO 
OO 2 N=1,NN 

CVX=0.0 

CYVY—On0 

CVZEO © 

NTEM=NN-N+ 1 

OO 3 LL=1,NTEM 
Li=LL+1 


10 


GAME VGRSEG(AF,YF.ZF.XTIPC(N,LL).YTIP(N,LL),ZTIP(N.LL) 
Seema ) ¥TLP(N, U1), ZTIP(N, L1), VX,VY.VZ.S%.SY,.SZ.0) 


CVX =CVX+VX 
CVY*CVY+VY 


3 CVZ=ECVZ+VZ 


VIX=VIX+GB(N,MM,K)*CVX Gi y 
VIY=VIY+GB(N,MM,K)#CVY 


2 VIZ=VIZ+GB(N,MM,K)*CVZ 
----- TRANSITION WAKE--------------------------------------------------- 


WAKX=0Q. 
WAKY =O. 
WAKZ=0. 


00 35 N=1,100 


CALL VORSEG(XF,YF,ZF,XH(N,1),YH(N, 1), ZH(N, 1), XH(N41, 


qo) ZHCNS 4, 1), VX.VY.VZ.SX.SY.SZ.O) 
WAKX =WAKX+VX 
WAKY =WAKY+VY 
35 WAKZ=WAKZ+VZ 
CALL VORSEG(XF,YF,ZF,XH(1,1),0..0.,XH( 101, 
& HBWX,HBWY,HBWZ,.SX,.SY,SZ.0) 
HBWX=-HBWX 
HBWY =-HBWY 
HBWZ=-HBWZ 
41. OO 71 M=1,MM 
IF(M.EQ.1) GO TO 72 
OO 73 N=1,NW 
UTX(N, 1) 2UTX(N, 2) 
UTY(N, 1) 2UTY(N, 2) 
UnZzCN 1 =UTZ(N. 2) 
OO 74 N=1,NW 


190 O% « 


is 
16? 


GALE VORSEG(XEF. YF,ZF,XW(N,M),YWC(N,M), ZWON,M) ,XW(N, oo YW(N,M+1), 


1 ZW(N,M+1),UBX(N),UBY(N),UBZ(N),SX.SY,SZ.0) 
IF(M.GT.1) GO TO 75 


CALL VORSEG(XF,YF.,ZF.XW(N,M),YW(N,M), ZWOIN,M),XW(N+1,M),YWO(Nt+1,M), 


jew Ne tem) oT X(N 1) 2UTYON, 1) UTZUN. 1).SX.SY.SZ.O) 


1), YHON+1, 


We 


75 CALL VORSEG(XF,YF,ZF,XW(N,M+1),YWON,M4+1), ZW(N,M+1) , XWO(N+1,M41), 
1 YW(N+1,M4+1),ZWON+1,M4+1),UTX(N,2),UTY(N,2),.UTZ(N,2),SX,SY,SZ.0) 


74 CONTINUE 
DO 76 N=1,NW 


C #** COMPUTE INOUCED VELOCITY FROM ONE FORCESHOE AT THIS RAOIUS 


C 


BUG=UTX(N,2)-UTX(N, 1)+UBX(N) 
CAT2UTY(N,2)-UTY(N, 1)+UBY(N) 
OOG=UTZ(N,2)-UTZ(N, 1)+UBZ(N) 
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UEFPOlS49 
Op ea @) 2/5) S1@) 
UR FOSS? 
UFPO1S570 
UFPO1580 
UrPPOIS90 
UFPO1600 
UFPO1610 
UEP O1620 
UFPO1630 
UFPO1640 
UFPO1650 
UFPO1660 
UFPO1670 
UFPO1680 
UFPO16390 
UFPO17CO 
UPFOI? 10 
URE Ol 2o 
UFPO1730 
UFPO1740 
URROI7 50 
UFPO1760 
UrRPOl? 70 
UFPO 1780 
UFPC17S0 
UFPO1800 
UFPO1810 
UFPO1820 
UFPO1830 
UFPO1840 
UEPOTS¢ 
UFPO1a60 
UFPO1870 
UFPO1880 
URFOVEsO 
UFPO1900 
UFPOQ1910 
UrFrPO1s2° 
UFPO19S0 
UFPO1940 
UEPOTI20 
UFPOQ1960 
WEEOGIS7TO 
UFPO1980 
UFFPOISSO 
UFPO2000 
UFPO2010 
UFPO2020 
UFPOQ2030 
UFPO2040 
UE PO2ZOS0 
UFPO2060 
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Ug! 


78 


va 
2 


93 


71 
10 


IF(N.NE.NW) GO TO 77 
IF(M.NE.MID) GO TO 78 


ALL TRANSITION WAKE VORTICES DISAPPEAR EXCEPT MID, 


WHICH TURNS 


VEPOZG7SO 
UFPO2080 
UF PO2090 
UFPQ2 100 


INTO ULTIMATE WAKE VORTICES (TIP ANO HUB) BY ASSUMING MEAN LOADINGUFPO2110 


BUG =BUG+WAKX+HBWX 
CAT=CAT+WAKY+HBWY 

DOG=DOG+WAKZ+HBWZ 

GO TO 78 

BUG=BUG-UBX(N+1) 

CAT=CAT-UBY(N+1) 

DOG=DO0G-UBZ(N+1) 

VOFWX(N,M)=BUG 

VOFWY(N,M)=CAT 

VOFW2Z(N,M)=D0G 

CONTINUE 

DO 93 LN=1,NVFW 

KKM=K-LN 

IF (KKM.LE.0O) KKM=KKM+4+NSR 
VIX=VIX+VOFWX(LN,M)=(GMEAN(M) 

1 +OKFW(LN)*(GAM(M,KKM) -GMEAN(M) )) 
VIY=VIY+VOFWY(LN,M)#*(GMEAN(M) 

1 +DKFW(LN)=(GAM(M,KKM)-GMEAN(M) ) ) 
VIZ=VIZ+VOFWZ(LN,M)=(GMEAN(M) 

1 +DKFW(LN)*=(GAM(M,KKM) -GMEAN(M) ) ) 
CONTINUE : 

RETURN 

END 
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UFPO2140 
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UFPO2160 
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UFPO2210 
URPO2220 
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UFPO2260 
UFPO2Z270 
UF PO2280 
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UEPO2Z329 
UFPO2330 
UFPO2340 
UFPO2Z350 
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SUBROUTINE RDFILE 


Sere eeaeeeazneReeeeneneecneaeenReeeeanszseo ee eeeneeweeneeeeaznenzneueeeneenee 


* 2 
. SUBRDUTINE RDFILE READS DISK FILE14 DATA CREATED * 
* BY PUF-2 AND XAVES DATA NEEDED BY FPUV IN _ 
* COMMONS GEDCOM AND SINCOM. NOTE THAT NSR IN * 
ce PUR 2eMusn SE DIVISTBEE SY THE NUMBER DF SLADES, * 
* NBLADE. i 
e * 
= * 
NN => = 
* MM => x 
= NW => NUMBER OF VORTICES IN KEY WAKE * 
= NBLADE => NUMBER DF PROPELLER BLADES - 
am NVC 22 No DF CHORD VDRTICES, OTHER BLADES * 
* NVS => NO OF SPAN VORTICES, OTHER BLADES : 
= NVW => NO OF WAKE VORTICES, OTHER WAKE : 
~ NSR 2=> NUMBER DF TIME STEPS/REVOLUTION be 
* NWK => NUMBER DF RADII IN WAKE HARMONICS : 
2 2 
Serer Reese eeeeeesteeeeseteerenee ee eeseeseereveee zene dd 


COMMON TD CONTAIN GEOMETRIC PARAMETERS 


COMMON / GEOCDM / MM,MN,NN,NW,NTORL ,X(12,10),Y(12,10).2(12,.10), 
& Sarre OM) Yate Ol1O.11).z271e( 10, 171), XH0 101.7) .YHC101,7). 
& ZH(101,7),.XW(21,10),YW(21,10),ZW(21,10),RH 


COMMON TO CONTAIN SINGULARITY STRENGTHS 


COMMON / SINCOM / GT(10,10,60),GB8(10.9.60).GTW(10),S8(10,9),GTV 
& ,GAM(9,60) ,GMEAN(9),OKFW(20) ,HBWX,HBWY ,HBWZ 


COMMON TO CONTAIN FIELO POINT DATA, N&SLADE, IDENT, NSR, AND VSR 


COMMON / FPCOM / NSLADE,NSR,XFP,RFP,TZFP,YBAR, ZBAR, 
& DSLADE ,OTFP,IDENT(18),TFP(60),VSR,IPLOT.DIAM 


DIMENSION R2Z(10),R(9),CHDCP(S) ,CHORZ(10) ,PHI8(90) 

DIMENSION X0(6,4,6),Y0(6,4,6),2Z0(6,4,6) 

DIMENSION RWK(11),AWA(16,11),8WA(16,11),AWR(16,11),8wR( 16,11), 
& AWT(16,11),.8WT(16,11) 

DIMENSION XWO(S.4,6),YWO(5,4,6),ZWO(5,4,6) 

DIMENSION S$80(5,3),G68C(5,3,60) ,GAMWC(4,3,60),GTC(5,4,60), 
& GAMC(3,60),GMEANO(3) ,OKOW(4) | 

READ(14) AJ 

IF (AJ.GE.2.0) STOP 

READ(14) ISTDY,NBLADE,MM,NN,NW,NVC,NVS,NVW,NSR,NWK 


SET CONSTANT FOR GRID 


= 70 = 


UF PO2360 
UEre237©0 
UF PO2380 
UFP@2 320 
UF PO2400 
UFPO2410 
UFPOZS20 
UFPO2430 
UFPO2440 
UFPO2450 
UF PO2460 
UFPO2470 
UFPO2480 
UFPO2490 
UF PO2S500 
UEEOZo1© 
UF PO2520 
UFPO2530 
UFPO2540 
UFPO2550 
UFPO2560 
UEPO2Z570 
UFPO2580 
UFPO25390 
UF PO2600 
UFPO2610 
UF PO2620 
UFPO2630 
UFPO2640 
URFOZ650 
UFPO2660 
UFPO2670 
UFPO2680 
UF PO2690 
UFPO2700 
UFPQ2710 
WIR PO2Z720 
UFPO2730 
UFPO2740 
UF PO27 50 
UFPO2760 
WEP @2 070 
UFPO2780 
UFPO2790 
UF PO2800 
UFPO28 10 
UF PO2820 
UFPO2830 
UFPO2840 
UFPO2850 
UFFPO28560 
UFPOQ2870 
UFPO2830 
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aan 
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onan ra WN 
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zs 


& WN 


NM=NN- 4 
MN=MM+ 4 

NP =NN+ 14 
NZ=NW+ 1 
NVCP=NVC+1 
NVSP=NVS+1 
NBLO=NBLADE- 1 


SET CONSTANTS FOR UNSTEADY PROBLEM 


IF (ISTOY.EQ.1) NSR=1 
NVYP =NVW+ 1 
NTB=NN*MM 


READVKEY SEADE GEOMETRY 


READ(14) ((X(N,M),N=1.NP),M=1,MN), 
CCY(N,M),N=1,NP),M2=1,MN), 
((Z(N.M).N#1,.NP),M=1,MN), 

. (RZ(M),M=1,MN), 
(R(M),M=1,MM), 
(CHOCP(M),M=1,MM), 
(CHORZ(M),M=1,MN), 
(PHIB(K),K=1,NTB) 


READ BLADE GEOMETRY OF DTHER BLADES 


Tee NSEADE SEE i)GeTO 11 


READ(14) (((XO(N,M,K).N=1,NVCP),M=1,NVSP),K=1,NBLD), 
((CYO(N,M,K),N=1,NVCP),M=1,NVSP),K=1,NBLO), 
(((ZO(N,M,K),N=1,NVCP),M=1,NVSP),K=1,NBLD) 


CONT INUE 
READ WAKE RADII AND WAKE HARMONICS 


READ(14) (RWK(M),M=1,NWK) 

READ( 14) ((AWA(I,M),I=1,16),M=1,NWK), 
((BWA(I,M),I=1,16),M*1,NWK), 
(CAWR(I,M),I12=1,16),M=1,NWK), 
((B8WRCI.M).I=1,16),M21,.NWK), 
(CAWT(I,M),I=1,16),M=1,NWK), 
((BWT(I,M),1=1,16),M=1,NWK) 


READ VARIOUS PROPELLER/WAKE CONSTANTS 


READ(14) WAKE,RPM,OIAM,RH,CDRAG,VSR,UR, (IOENT(N) ,N=1, 18) 


READ KEY WAKE GEOMETRY 


((XW(N,M),N21,NZ),M=1,MN), 
(CYW(N,M),N=#1,NZ).M=1,MN), 
((ZW(N,M),N=1,NZ),M=1,MN), 
((XH(N,K),N#1, 101),K=1,NBLADE), 


READ( 14) 


a 


UFPO2890 
UFPO23900 


« UERO2ZS 10 


UPPO2Z92° 
UFPO2930 
UFPO2940 
UFPO23950 
UFPO2960 
UFPO2970 
UFPO2980 
URPO23az© 
UF PO3000 
UFPO3010 
UF PO2ZO2Z¢ 
UF PO3030 
UFPO3040 
UFPOSOS0 
UFPO3060 
UFPO3070 
UFPO3080 
UFPe2C39o 
UFPO3 100 
UFPO3110 
UFPOS 120 
UFPO3130 
UFPO3 140 
UFPO3 150 


‘UFPOS160 


UFPO3170 
UFPO3180 
UFPO3190 
UFPO3200 
UFPO3210 
UFPO3220 
UFPO3230 
UFPO3240 
UFPO3250 
UFPO3260 
UFPO3270 
UFPO3280 
UFPO3230 
UFPO3300 
UFPO3310 
UF PO3320 
UFPO3330 
UF P03340 
UFPO3350 
UFPO3360 
UFPO3370 
UF PO3380 
UFPO3390 
UFPO3400 
UFPO3410 





Cc 
C 
C 


Cc 
Cc 
Cc 


Cc 
C 
C 


C &&& 


((YHCN,K),N=1,101),K=1,NBLADE), 
((ZH(N,K).N=1,101),K=2=1,N8LADE), 
((XTIP(N,L),L21,NP),N=1. NN), 
(CYTIP(N,L),L=1,NP),N=1,NN), 
((ZTIP(N,L),LA1,NP),N=1,NN), 


omon au 


sae READ KEY BLADE SINGULARITY STRENGTHS 


((SB(N,M),.N=1,.NN),M=1,MM), 
(((GB(N,M,KK),N=1,NN),M=1,MM),KK=1,NSR), 
(( (GT(N,M,KK),N=1,NN),M=1,MN),KK=1,NSR), 
((GAM(M,KK),M=1,MM),KK=1,NSR) 


WN — & 


wae OTHER WAKE GEOMETRY 


IF (NBLADE.LE.1) GOTO 12 

READ(14) (((XWO(N.M,K).N=1,NVWP),M=1,NVSP),K=1,NBLO), 
2 ((CYWO(N,M,K),N=1,NVWP),M=1,NVSP),K=1,NBLO), 
3 (((ZWO(N.M,K).N=1,NVWP),M=1,NVSP),K=1,NB8LO), 


xee2 OTHER BLADE SINGULARITIES 


2 
IS) 


901 
Oz 
903 


((SBO(N,M),N=1,NVC),M=1,NVS), 
(((GBC(N,M,KK),N=1,NVC),M=1,NVS),KK=1,NSR), 
(( (GAMWC(N,M, KK) ,N=1,NVW),M=1,NVS),KK=1,NSR), 
(CC GTCC(N,M,KK),.N=1,NVC),M=1,NVSP),KK=1,NSR), 
((GAMC(M,KK),M=1,NVS),KK=1,NSR), 
(GMEAN(M),M=1,MM), 

(GMEANO(M) ,M=1,NVS), 

(DKFEFW(N),N=1,NW), 

(DKOW(N) ,N=1,NVW) 


NYrwaaewM Orang iuw»nh 


CONT INUE 
GOTO 100 


DIAGNOSTICS &&&&&8384444834558538488& 


FORMAT( 1615) 
FORMAT(8F 10.5) 
FORMAT(18A4) 
WRITE( 1,902) AJ 
WRITE( 1,901) ISTDY,NBLADE,MM,NN,NW,NVC,NVS,NVW,NSR, NWK 
WRITE( 1,902) ((X(N.M),N=1,NP),M=1,MN), 
CCYO(N,M) .N2=1,NP),M=1,MN), 
((Z(N,M),N21,NP),M=1,MN), 
(RZ(M),M=1,MN), 
(R(M),M=1,MM), 
(CHOCP(M),M=1,MM), 
(CHORZ(M),M=1,MN), 
(PHIB(K),K21,NTB) 
IF (NBLADE.LE.1) GOTO 21 
WRITE( 1,902) (((XO(N,M,K),.N=2=141,NVCP),M=1,NVSP),K=1,NBLO), 
2 ((CYO(N,M,K),N=1,NVCP),M=1,NVSP),.K=1,NBLO), 
3 (CC ZO(N,M,K),N21,NVCP),M=1,NVSP),K=1,NBLO) 


Ono uh Wh 


DE 


UFPO3420 
UFPO3430 
UFPO3440 
UFPO3450 
UFPO3460 
UFPO3470 
UFPO3480 
UFPO3490 
UEP O3500 
UFPOI5S 10 
UEPO2520 
UPPOsSS20 
UFPO3540 
UF POSS350 
UF POsS00 
UFPO3570 
UFPO3580 
UFPO2S9S0 
UFPO3600 
UFPO36 10 
UFPO3620 
UFPO3630 
UFPO3640 
UVEFOS630 
UFPO3660 
UFPOS6E70O 
UFPO3680 
VErPOse sO 
UFPO3700 
UFPO37 10 
UF PO3720 
UFPO3730 
UEP O3740 
UFPO3750 
UFPO3760 
UFPO3770 
UFPO3780 
UFPO3790 
UFPO3800 
UFPO3810 
UFPOSE2Z0 
UFPO3830 
UFPO3840 
UFPO3850 
UFPO3860 
UFPO3870 
UFPO3880 
UFPO3890 
UF PO3900 
UFPO3S9I10 
UFPO3920 
UFPO3930 
UFPO3940 





21 CONTINUE UErPO3IS¢0 


WRITE( 1,902) (RWK(M) ,M=1,NWK) UFPO3960 
WRITE( 1,902) ((AWA(I.M),1=1,16),.M=1,NWK), UF PO3970 
3 ((BWA(I,M).1=1,16),M=1,NWK), UF PO3380 
4 ' ((AWR(I,M),1=*1,16),.M=1,NWK), UFPO3990 
5 ((BWR(I.M).1=1,16),M=1,NWK), UFPO4000 
6 ((AWT(CI,M),1=1,16).,M=1,NWK), UFPO4010 
7 ((BwWT(I,M),1=1,16).M=1,NWK) UF PO4020 
WRITE( 1,902) WAKE,RPM,DIAM,RH,CORAG,VSR,UR UFPO4030 
WRITE(1,903) (IDENT(N) ,N=1, 18) UFPO4040 
WRITE( 1,902) ((XW(N,M) .N=1,NZ),M=1.MN), UFPO4050 
2 ((YW(N,M) .N=1,NZ),M=1,MN), UFPO4060 
3 (( ZW(N,M),N=1,NZ),M=1,MN), UFPO4070 
4 ((XH(N,K),N=1,101),K=1,NBLADE), UFPO4080 
5 (CYH(N,K),N#1,101),K=1,NBLADE), UFPO4090 
6 ((ZH(N,K).N=1,101),.K=1,NBLADE), UFPO4 100 
ys ((XTIP(N,L),L=1,NP).N=1,NN), UFPO4110 
8 (CYTIP(N,L),L=A1,NP),N=1,NN), UFPO4120 
9 ((ZTIP(N,L),L=1,NP),N=1.,NN) UFPO4130 
WRITE( 1,902) ((SB(N,M).N=1.NN),M=1,MM), UFPO4140 
2 (((GB(N,M,KK) .N=1,NN),M=1,MM),KK=1,NSR), UFPO4150 
3 (((GT(N,M,KK),N=1,NN),M=1,MN) ,KK=1,NSR), UFPO4 160 
4 ((GAM(M,KK),M=1,MM),KK=1,NSR) UFPO4170 
IF. (NBLADE.LCE.1) GOTO 22 UFPO4180 
WRITE( 1,902) (((XWO(N,M,K),N=1,NVWP),M=1,NVSP),K=1,NBLO), UFPOQ4190 
D, (( CYWO(N,M,K).N=1,NVWP) .M=1,.NVSP),.K=1,N8LD), UF PO4200 
3 (((ZWO(N,M,K),N=1,NVWP),M=1,NVSP),K=1,NBLO) UFPO4210 
WRITE( 1,902) ((SBO(N.M) ,N=1,NVC),M=1.NVS), UFPO4220 
2 (((GBC(N,M,KK),N=1,NVC),M=1,NVS),KK=1,NSR), UF PO4230 
3 (((GAMWC(N,M,KK),.N=1,NVW) .M=1,.NVS),.KK=1,NSR), UFPO04240 
4 (((GTC(N,M,KK) ,N=1 NVC), M=1 -NVSP), KK=1, NSR), UFPO4250 
S ((GAMC(M,KK) ,M=1 -NVS). KK=1 _NSR)., UFPO4260 
5 (GMEAN(M) ,M=1, MM), UFPO4270 
7 (GMEANO(M),.M=1,NVS), UFPO4280 
8 (DKFW(N) ,N=1-NW), . UFPO4290 
9 (DOKOW(N) .N=1,NVW) UFPO4300 
22 CONTINUE UFPO4310 
100 RETURN UFPO4320 
END UFPO4330 
Foe 


-73- 
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SUBROUTINE ROTERM(ISTOP) 


Reese aanaueeee auc Bau egegngge ae eee eee eee eee ee ee ee Ree Vee eee Se 


= 
a 
x 
a 
* 
a 
= 
a 
”" 
*“ 
x 
x 
* 
x 
a 
Zz 
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SUBROUTINE ROTERM REAOS IN BLAOE ANO GEOMETRY 
OATA INTERACTIVELY ANO COMPUTES NECESSARY FIELO 
POINT GEOMETRY. 


NBLAOE 
NTORL 


XFP 


RFP 

TZFP 
YBAR, ZBAR 
TH 

Y8,Z8 


= 0 
= 1 
2 2 


NUMBER OF BLAOES 

CONSIOER THICKNESS ANO LOAOING 
CONSIOER THICKNESS ONLY 
CONSIOER LOADING ONLY 

AXIAL FIELO POINT COOROINATE, 
POSITIVE OOWNSTREAM 

RAOIAL FIELD POINT COORDINATE 
INITAL FIELO POINT ANGLE IN OEG. 
Y ANO Z FIELD POINT COOROINATES 
ANGLE FROM TIME STEP INPUT 

Woo leeOs PT eCOORO, SSe50RE TIME S 


Bueurtruaeeuaceeegeeueneueeeneneve eee eeezgceeveeese eee neeseee eee Reeesungeae see 


COMMON TO CONTAIN GEOMETRIC PARAMETERS 


T 


¢* #¢ #¢ # © © # @ © © @ @ @ &@ © & 8 8 


COMMON / GEOCOM / MM,MN NN ,NW,NTORE,X(12,.10),7(12, 10),2(12, 10), 


COMMON TO CONTAIN FIELO POINT OATA, 


NBLAOE, IOENT, 


COMMON / FPCOM / NBLAOE,NSR,XFP,RFP,TZFP,YBAR,ZBAR, 


OBLAOE ,OTFP,IOENT(18),TFP(60),VSR,IPLOT,OIAM 
FORMAT(/’ ENTER XFP, RFP, TZFP.// 
FOR HELP, ENTER -11. FOR XFP.°’) 


WRITE(6,901) 


READ(S,*«) XFP,RFP,TZFP 
IF (XFP.LT.-10.) GOTO 30 
{F (TZEP-LT 20-0) GOTO 70 
IF (RFP.GT.O.) GOTO 60 
WRITE(6,902) 


FORMAT(/? 


ENTER NBWLN,M,KK’/% 


FOR MECr, 


REAO(S,*) NBW,N,M,KK 


TH=6.283185*FLOAT(KK-1)/FLOAT(NSR) 


IF (NBW.EQ.0) GOTO 80 


FIELO POINT CENTEREO ON 


BLAOE GRIO. 


IF (NBW.E£Q.1) CALL FPGRIO(X,Y,Z,N.M,KK, 12) 


FIELD POINT CENTERED ON TRANSITION WAKE GRIO. 


IF(NBW.EQ.2) CALL FPGRIO(XW,YW,ZW,N,M, KK, 21) 


Say ie 


NSR, 


TPueelort terre (lO )) 2hlPpt 10.411) .XH( 101,7), YH( 101.7). 
ZH(101,7),XW(21,10),YW(21,10),2ZW(21,10),RH 


ANO VSR 


ENTER ZERO FOR NBW.‘) 


UFPO4340 
UFPO4350 
UFPO4360 
UFPO4370 
UFPO4380 
UFPO43390 
UFPO4400 
UFPO4410 
UFPO4420 
UFPO4430 
UFPO4440 
UFPO4450 
UFPO4460 
UFPO4470 
UFPO4480 
UFPO4490 
UFPO4500 
UFPO4510 
UFPO4520 
UFPO4530 
UFPO4540 
UFPO4550 
UFPO4560 
UFPO4570 
UrPO04580 
UFPO4590 
UF PO4600 
UFPO4610 
UFPO4620 
UFPO4630 
UFPO4640 
UFPO4650 
UFPO4660 
UFPO4670 
UFPO4680 
UFPO4630 
UFPO4700 
UFPO4710 
UFPO4720 
UFPO4730 
UFPO4740 
UFPO4750 
UFPO4760 
UFPO4770 
UFPO4780 
UFPO4730 
UFPO4800 
UFPO4810 
UFPO4820 
UFPO4830 
UFPO4840 
UFPO4850 
UFPO4860 





60 
70 


80 
903 


30 
304 


é 


& 
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& 
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& 
& 
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IF(NBW.NE.1.ANO.NBW.NE.2) GOTO 10 UF PO4870 
OTFP=360.0/FLOAT(NSR) UFPO4880 
RETURN UFPO4890 
ISTOP=1 UFPO4900 
RETURN UFPO4910 
WRITE(6,9303) UFPO4920 
FORMAT(/’ NBW = O => HELP’/S5X,’= 1 => BLADE WAKE GRID’/S5X, UFPO4930 
‘= 2 => TRANSITION WAKE GRIO’/’ N’,7X,%=> CHOROWISE COOROINATE’/UFPO4940 
’ M’ .7X,’=2> SPANWISE COORDINATE’/’ KK’,6X, UFPO4950 
‘=> TIME STEPS ROTATEO + 1’/) UFPO4960 
GOTO 20 i UFPO4970 
WRITE(6,904) UFPO49380 
FORMAT(/’ XFP <-10=> HELP’/8x, UFPO4390 
’>-10=> AXIAL FIELO POINT COOROINATE, POSITIVE OOWNSTREAM’/ UFPOSOOO 
‘ REP < © => SELECT COORDINATE USING BLADE OR WAKE GRIO’/ UFPOSO10 
8X,’> O => RADIAL FIELD POINT COORDINATE’/ UFPOSO20 
’ TZFP < 0 => STOP’/8X,’> O => ANGULAR FILE POINT COORDINATE’ UFPOSO30 
/) UFPOSO40 
GOTO 10 UFPOSOSO 
ENO UF POSO6O 


=) 7 
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SUBROUTINE FPGRIO(XX,YY.ZZ.N,M,KK,NDOIM) 


B@Beeeeneeeneeaeeene ee eee eae eeeweaeeee zee eee eeee ee eezexne eevee ensaeewze ee 8 
x 
SUBRDUTINE GRID CALCULATES X,R, AND THETA FIELD 
POINT GEOMETRY COMPONENTS GIVEN GRIO VECTORS 
X.Y, ANO Z FOR BLAOE OR WAKE, GRIOQ CDDROINATES 


N AND M, AND TIME STEP KK. 


* 

bg 
* @ 
* ” 
= 2 
x 2 
* XXA.NY.22 22> GRIO GEDMETRY VECTORS = 
* N,M => GRIOINOIGES SEEECTED FOR VELDCITY ~* 
* KK => TIME STEP. NUMBER OF ANGLE STEPS + 1* 
* XFP,Y8AR, ZBAR,TZFP,RFP * 
* pox Y ve crrETA.RAOIAL FIELD PT coun « 
* => ANGLE BETWEEN BLAOES * 
* => ANGLE OF TIME STEP 2 
* => FIRST OIMENSION DF GRIDS PASSED * 
2 = 
x 2 


OBLADE 
10: 
NOIM 


RRS RETR RERROKRHReAaERAaR ee eee eevee eeaeegeeeaeezge eee eae eee ae eee ee ee 


COMMON TO CONTAIN FIELD PDINT DATA, NBLAGE, IDENT, NSR, 

COMMON / FPCDM / NBLADE,NSR,XFP,RFEP,TZFP, YBAR, ZBAR, 
OBLAOE,OTFP,IDENT(18),TFP(60),VSR,IPLOT,DOIAM 

OIMENSION XX(NDIM,10),YY(NDIM, 10), ZZ(NOIM, 10) 


COMPUTE X,Y,Z CDDROINATES IN ZERO ANGLE POSITION 
TH=6.283185*FLDAT(KK-1)/FLDAT(NSR) 
THOEG=TH*360./6.283185 

OBLAQE=360.0/FLDAT(NBLAOE ) 

XFP=0.25*(XX(N,M) + XX(N#1,M) + XX(N,Mt+1) + XX(N4+14,M41)) 
YB =0.25*(YY(N,M) + YY(N+1.M) + YY(N,M41) + YY(N+4,M+1)) 
ZB =O.25*(ZZ(N,M) + ZZ(N+1,M) + ZZ(N,M4+1) + ZZ(N+1,M+1)) 


CONVERT TD CYLINDRICAL CDDRDINATES 


RFP=SORT(Y8*Y8+Z8#2Z8) 
TZFP=ATAN2(Z8,Y8)/1.7453293E-02 
TZEPSTZEP+TRHOEG 


BRANCH IF GENERATING PIEWAKE OATA 


IF (IPLOT.EQ.2) GOTO 200 


TZFP=180.0/FLDAT(NBLAQE) - ATAN2(Z8,YB)/1.7453293F-02 


KEEP TZFP BETWEEN O AND DONE BLAOQE ANGLE 
OD 100 K=1,NBLAOE 

IF(TZFP.LT.0.0) TZFP=TZFP+OBLAOE 
IF(TZFP.GT.OBLAOE) TZFP=TZFP-DBLAOE 
CONTINUE 
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ANOQ VSR 


UFPOQSO70 
UF POSO80O 
UFPOSO3S0 
UFPOS 100 
UFPOST10 
UFPOS120 
UFPOS130 
UFPOS140 
MER@S 130 
UFPOS160 


~UFPOS 1702 


UF POS 180 
UFPOS190 
UFPO5200 
UFPOS210 
UFPOQS220 
UFPOS230 
UFPOQS240 
VEPGs25¢0 
UFPOQS260 
UPPOS27¢ 
UF POS280 
UFPOS290 
UFPOS300 
UFPOS53 10 
UF POS320 
UFPOS330 
UFPOQS340 
UF POSS3S50 
UFPOS360 
UFPOS370 
UFPOS380 
UFPOS390 
UFPOS400 
UFPOQS410 
UFPOS420 
UFPOQS430 
UFPOS440 
UFPOQS450 
UF POS460 
UFPOS470 
UFPOS480 
UFPOS490 
UF POSS00 
UFPOSS10 
UPPGSs2¢0 
UEPOGSs2C 
UFPOSS40 
UFPOSSSO 
UFPOSS60 
UFPOSS7¢0 
UFPOS580 
UFPOSS390 





C sue 


200 


300 


ROTATE BLADE FOR PIEWAKE DATA 


YBAR=YB*COS(TH) - Z8*SIN(TH) 
ZBAR=YB*SIN(TH) + ZB*COS(TH) 
RETURN 

END 


ee ae 


UFPO5S600 
UFPOS610 
UFPOS620 
UFPOS630 
UFPOS5640 
UFPOS5650 
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SUBROUTINE FPSTEP(K) 


aanraneka eee eeaeaan Rea eaeeaeagneeaeageeaeaeaeneeeaeeaeaeenaecneaeeeceaczne eee ee 


x x 
* FPSTEP IS A SUBROUTINE TO CALCULATE THE FIELO * 
* POINT VELOCITIES INOUCEO BY ONE BLAOE AT “TIME * 
* STEP*, K. THE "“OIRTY WORK" IS OONE IN * 
* SUBROUTINE FPVEL. * 
= = 
* TFP => ANGULAR FIELD POINT COOROINATE, OEG. * 
* TRAO => ANGULAR FIELO POINT COOROINATE, RAO. * 
* YFP,ZFP => Y ANO Z FIELO POINT COOROINATES * 
* VIX,VIY,VIZ,VIR,VIT => INOUCEO VELOCITY COMPONENT = 
x * 
eee aeaeecznaczgnea Rae eee zeeeeazgnaeneeaeeeaeeeczneeeaeaeaneanaenennnnnaE 


**= COMMON TO CONTAIN GEOMETRIC PARAMETERS 

COMMON / GEOCOM / MM,MN,NN,NW,NTORL,X(12,10),Y(12,10),2(12,10), 
& ATIPC1O- t1) VY TIPt tO, (1), ZT1P (10.11). XH( 101.7). YHC101,7). 
& ZH(101,7),XW(21,10),YwW(21,10),2W(21,10)},RH 


sxe COMMON TO CONTAIN SINGULARITY STRENGTHS 


COMMON / SINCOM / GT(10, 10,60),G8(10,9,60) ,GTW(10),S8(10,9),GTV 
&  ,GAM(9,60),GMEAN(9),OKFW(20),HSWX ,HBWY,HBWZ 
**= COMMON TO CONTAIN FIELD POINT OATA, NBLAOE, IOENT, NSR, ANO VSR 


COMMON / FPCOM / NBLACE,NSR,XFP,RFEP,TZFP,YBAR, ZBAR, 
& OBLAOE ,OTFP, IOENT(18),TFP(60),VSR,IPLOT,OIAM 


eae COMMON TO CONTAIN VELOCITY OATA 
COMMON / VELCOM / VIX(60).VIT(60),VIR(60),U(60,.3),NSAMP 
CALCULATE CARTESIAN FIELO POINT COOROINATES TO FEEO FPVEL 


TFP(K)=FLOAT(K-1)*OTFP-TZFP 
TRAD=TFP(K)*1.7453293E-02 
COSFP=COS( TRAD) 
SINFP=SIN(TRAO) 
YFP=RFP*COSFP 

ZFP=RFP*SINFP 


ANO VIZ FOR FIELO POINT 


wee FPVEL WILL CALCULATE VIX, VIY, 
GALL FREVEL(K.XFP,YFP.ZFE.VIX(K),VIY,.VIZ) 
VIR(K)=SVIY*COSFP + VIZ*SINFP 
VIT(K)=VIZ*COSFP - VIY*SINFP 

RETURN 

END 
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UF POS660 
UFPOS670 
UF POS680O 
UVEPOS6 90 
UFPOS700 
UFPOS7 10 
UFPOS720 
UFPOS730 
UFPOS740 
UFPOS750 


' UFPOS760 


UFPOS770 
UFPOS780 
UFPOS790 
UF POS800 
UFPOS810 
UFPOS8&20 
UFPOS830 
UFPOS840 
UFPOS8S0O 
UFPOS860 
UFPOS870 
UFPOS880 
UFPOS890 
UFPOSSOO 
UFPOS910 
UFPO5920 
UFPOS930 
UFPOS940 
UFPOS950 
UFPOS960 
UFPOS970 
UFPOS980 
UFPOS5390 
UF PO6O00O 
UFPO6010 
UFPO6020 
UFPO6030 
UFPO6040 
UFPO6O0SO 
UFPO6060 
UFPO6070 
UFPO6O080 
UFPO6090 
UFPO6 100 
UFPO6110 
UFPO6 120 
UFPO6130 
UFPO6 140 
UFPO6150 
UFPO6160 
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SUBROUTINE SUMVEL 


eseeeenaeeenee eae eee eaeneeeeaneaeaeeeeeeenegnene eee eee nee aczeeeannenes 


a 


* BLADE, 

* TOTAL 

* VIX,VIR,VIT => 
* 

*  U(N,J) => 
* 

* J=O => 
* J=1 => 
* J=2 => 
od NS AMP => 
* 

PY 

» 

* 


COMMON TO CONTAIN FIELO POINT DATA, 


COMMON / FPCOM / NBLAOE,NSR,XFP,RFP,TZFP,YBAR, ZBAR, 


SUBROUTINE SUMVEL SUMS VELOCITIES FROM EACH 
CALCULATED IN FRSTEF, 
PROPELLER-INDUCED FIELO POINT VELOCITY. 


TO” PRODUCE THE 


VELOCITY COMPONENTS INOUCED BY 
ONE BLAOE AT "TIME STEP” K 
VELOCITY INOUCED BY ALL BLAOES 
AT “TIME STEP" N 

AXIAL COMPONENT 

RAOTAL COMPONENT 

TANGENTIAL COMPONENT 

NUMSER OF "TIME STEPS” IN ONE 
BLAOE ANGLE (= NUMBER OF 

U VELOCITIES COMPUTED. ) 


a2eegQgeaetree eee eeenegReReeaBeeeean an eseeaeezae eeezgeegeeezeegeaeeg Ree aR eee ESE 


NBLAOE, 


OBLAOE,OTFP,IDENT(18),TFP(60),VSR,IPLOT,OIAM 


#** COMMON TO CONTAIN VELOCITY DATA 


“*"e 


10 


20 


wee 


30 


aé 


COMMON / VELCOM / VIX(60),VIT(60),.VIR(60),U(60,3),NSAMP 


AT ANY GIVEN REAL TIME, 
TRAVELLEO THROUGH A CERTAIN NUMBER OF 
BY SUMMING VIX, VIR, 


ANO VIT AT BLAOE RATE, 


EDEN, 


* 
bs 
* 
= 
* 
a 
* 
= 
* 
zm 
2 
= 
a 
ad 
" 
* 
* 


NSR, 


THE VECCCITIES INDUCED BY EACH BLAOE INTO ONE VELOCITY 
INOUCED BY THE ENTIRE PROPELLER. 


NS AMP =NSR/NBLAOE 
DO 20 N=1,NSAMP 
00 10 J=1,3 
U(N,J)=0.0 
DO 20 K=1,N8LAOE 
L=N+(K-1)*NSAMP 


U(N, 1) =U(N, 1 )4+VIX(L) 
U(N,2)=U(N,2)+VIR(L) 
U(N,3)=U(N,3)4#VIT(L) 


NOW CONVERT FROM REFERENCE VELOCITY TO SHIP VELOCITY 


00 30 N=1,NSAMP 
00 30 vF1,3 
U(N,J)sU(N,JU)/VSR 
RETURN 

ENO 


=- 79 - 


ANO VSR 


EACH BLAOE CAN BE CONSIDERED TO HAVE 
PRIME STcrs” . 
WE CAN INCLUOE 


UFPO6170 
UFPO6 180 
UFPO6190 
UFPO6200 
UFPO6210 
UFPO6220 
UFPO6230 
UFPO6240 
UFPO6250 
UFPO6260 
UFPO6270 
UFPO6280 
UFPO6290 
UFPO6300 
UFPO6310 
UFPO6320 
UFPO6330 
UFPO6340 
UFPO6350 
UFPO6360 
UFPO6370 
UFPO6380 
UFPO6390 
UF PO6400 
UFPO6410 
UFPO6420 
UFPO6430 
UFPO6440 


UFPO6450 


UFPO6460 
UFPO6470 
UFPO6480 
UFPO64390 
UFPO6500 
UPPOSS 10 
UFPO6S520 
UFPO6S30 
UFPO6540 
UFPO6S50 
UFPO6560 
UFPO6S7/70 
UFPO6580 
UFPOGSSO 
UF PO&66EOO 
UFPO66 10 
UFPO6620 
UFPO6630 
UFPO6640 
UFPO66E50 
UFPO6660 
UFPO6670 
UFPO6680 
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SUBROUTINE FOUCAL UFPO6690 


RERRREEREREHRRERSER ERR ERE REM eee eee eee keke eee eek UFPO67O00 
* = UFPOQ6710 
* SUBROUTINE FOUCAL TAKES THE PROPELLER INOUCED * UFPO6720 
* VELOCITIES CALCULATEO IN SUMVEL AND PERFORMS * UFPO6730 
* A HARMONIC ANALYSIS ON THEM. IT CALCULATES * UF PO6740 
* FOURIER COEFFICIENTS FOR BOTH COSINE (A) ANO = UFPO67S50 
* SINE (8) ANO ALSO COSINE AMPLITUOE (AMP) AND * UFPO6760 
* PHASE ANGLE (PH). * UFPO6770 
* * UFPO6780 
* NH => NUMBER OF HARMONICS CALCULATED * UF PO6790 
* A(N,J).8(N,JU) => FOURIER COEFFICIENTS a UFPO6800 
* A => COSINE COEFFICIENTS 2 UFPO6810 
* B => SINE COEFFICIENTS * UFPO6820 
* N => HARMONIC NUMBER UP TO 15 * UFPO6830 
* vez => AXIAL VELOCITY * UFPO6840 
* => 2 => RAOIAL VELOCITY * UFPO6850 
= = 3 => TANGENTIAL VELOCITY * UFPO6860 
* AMP(N,J) => FOURIER COEFFICIENT AMPLITUDE * UFPO6870 
* N AND J SUBSCRIPTS AS FOR A, B ®* UF PO6880 
* PH(N,J) => PHASE ANGLE ASSOCIATEO WITH AMP * UFPO6890 
* NFB => NUMBER OF FIELO POINTS / BLAOE * UFPO6900 
* * UFPO6910 
Regge eeeaunueeceeeune eueeeeeeeueaueaneeueeeueneueeeeeeeneeeeeeaneeunee UF PO6920 
UF PO6930 

*#* COMMON TO CONTAIN FIELO POINT OATA, NBLADE, IDENT, NSR, ANO VSR UF PO63940 
UFPO6950 

COMMON / FPCOM / NBLAOE,NSR,XFP,RFP,TZFP, YBAR, ZBAR, UF PO6960 
& OBLAOE ,OTFP, IOENT(18),TFP(60),VSR,IPLOT,OIAM UFPO6970 
UFPO6980 

*** COMMON TO CONTAIN VELOCITY DATA UFPO63990 
UFPO7000 

COMMON / VELCOM / VIX(60),VIT(60),VIR( 60) ,.U(60,3),NSAMP UFPO7010 
: UFPO7020 

*a* COMMON TO CONTAIN HARMONICS COEFFICIENTS UFPO7030 
UFPO7040 

COMMON / HARCOM / A(15,3),8(15,3),AMP(15,3),PH(15,3),NH UFPO70S0 
UFPO7060 

eae FIRST CALCULATE A ANDO B UFPO7070 
- UFPO7080 

NFB =NSR/NBLAOE UFPO7090 
NHMAX=NFB/2+1 UFPO7 100 
NH=MINO(NHMAX, 15) UFPO7110 
STEP=6.28318S/NFB UFPO7120 
ANFB=2.0/NFB UFPO7130 
BNFB=1.0/NFB UFPO7140 
00 11 K=1,NH UFPO7150 
00 10 J=1,3 UFPO7160 
A(K,J)=0.0 UFPO7170 
10 B(K,J)=0.0 UFPO7180 
DO 9 N=1,NFB UFPO7190 
ee CUN= 1) =(K=1))*STEP UFPOQ7200 
ST=SIN(T) UFPO7210 


~ &0 « 





C «a 


a 


ue] 


11 


aé 


GT =COS (7 ) 


00 9 Jz1,3 
A(K,J)=A(K,J)+U(N, JU) *CT 


B(K,J)=B(K,J)tU(N,J)*ST 


CONVERT A AND B INTO AMP AND PH 


DO 11 J#1,3 


IF(K.EQ.1) AC 1,U)=A(1,JU)*BNFB 


IF(K.EQ.NHMAX) GO TO 


Ue 


IF(K.GT.1) A(K,J)=A(K,U) *ANFB 
IF(K.GT.1) B(K,J)=B(K,J)*ANFB 


GO TO 13 


A(K,J)=A(K,JU)*8NFB 

B(K,J)=0.0 
AMP(K,J)=SORT(A(K, J) **2+8(K,J) #22) 
PH(K,J)=ATAN2(A(K,JU),8(K,J))*57.296 
IF(PH(K,JU).LT.0.0) PH(K,J)=PH(K,JU)+360.0 


PH(K,J)=PH(K, J) /NBLADE 


PH(1,JU)20.0 
AMP(1,JU)=A(1,U) 


CONTINUE 
RETURN 
END 


Sacks 


UFPO7 220 
UFPO7230 
UFPO7240 
UFPO7250 
UFPO7260 
UFPO7270 
UFEPO72Z230 
UFPO729S0 
UFPO7300 
UFPO73 10 
UFPO7320 
UF PO7330 
UFPO7340 
UFPO7 350 
UFPO7360 
UFPO7370 
UFPO7380 
UFPO?7 330 
UFPO7400 
UFPO7410 
UFPO7420 
UFPO7430 
UFPO7440 
URFrPO7435°0 
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SUBROUTINE PIEDAT(ISTOP) 


aaeeneeeueeseeeuneazaeenaeeeeneneeencneaenzneeeznneaeeannee sn eecneenence eee 
= 

SUBROUTINE PIEOAT GENERATES OATA FROM INPUT TO 
AXIAL VELOCITY OUTPUT TO BE USEO AS INPUT FOR 
UNSTEAOY PIEWAKE. IT USES BLADE GRIO GEOMETRY 

TO CALCULATE NSR FIELO POINT AXIAL TIME-AVERAGED 
VELOCITIES AT FIVE RADII ALONG THE LEADING EOGE. 
IT THEN CALCULATES NSR FIELO POINT AXIAL TIME- 
AVERAGED VELOCITIES AT 1.1 TIMES PROPELLER RAOIUS 


N => CHOROWISE GRID COORDINATE 

M => SPANWISE GRIO COORDINATE 

K,KK => ANGLE INCREMENTS + 1 

ISTOP=1 => STOP FPUV UPON RETURN TO MAIN PROG 
NSR => NUMBER TIME STEPS IN REVOLUTION 


NR => RADIUS NUMBER 
XFP,YBAR,ZBAR,RFP,TZFP 
=> X,Y,Z,RAOIAL, THETA FIELO PT COORD 
THETA(1-5)=> ANGLES 5S TIME STEPS BEFORE ZERO 
BLAOE ANGLE 
THETA(6 - NSR+5) => ANGLES IN ONE BLAOE REV. 
THETA(NSR+6 - NSR+20) => ANGLES 15 TIME STEPS 
AFTER ZERO BLAO ANGLE. 


VX(K) => AXIAL VELOCITY AT ANGLE THETA(K) 
UA(KK,M) => SMOOTH AXIAL VELOCITY AT RADIUS M, 
ANGLE KK 


* £ *£ # © &# © &@ © © @ © © © @ © &@ © © © © © & 08 @ & 
¢ *#* @# &# © © © @ © © © © © © © #@ © © &@ © &@ @ @ 08 8 8 & 


zaeaxeeeaeenczneeeeeeeaaeeeeeeeneneeeneeneneneeeeeeneeneceeenene eee 


sae COMMON TO CONTAIN GEOMETRIC PARAMETERS 


¢ 


COMMON / GEOCOM / MM,MN,NN,NW,NTORL,X(12,10),.¥(12,10),2(12,10), 


& Pere ClO. yi lPt1O,11),Z251P (10, 11), XH( 101,7),YH( 101, 


& BO tet awial. 10), 1W( 21,10), 2Ww(21,10).RH 


*** COMMON TO CONTAIN SINGULARITY STRENGTHS 


COMMON / SINCOM / GT(10,10,60),G8(10,9,60),GTW(10),$8(10,9),GTV 


&  ,GAM(9,60),GMEAN(9),OKFW(20),HBWX,HBWY ,HBWZ 
**= COMMON TO CONTAIN FIELO POINT OATA, NBLADE, IOENT, NSR, 


COMMON / FPCOM / NBLAOE,NSR,XFP,REP,TZFP,YBAR, ZBAR, 
& OBLADE ,OTFP,IOENT(18),TFP(60),.VSR,IPLOT,DOIAM 


**« COMMON TO CONTAIN VELOCITY OATA 
COMMON / VELCOM / VIX(60),VIT(60),VIR(60),U(60,3),NSAMP 


wee COMMON TO CONTAIN HARMONICS COEFFICIENTS 


SED o 


We 


ANO VSR 


UFPO7460 
UFPO7470 
UFPO7480 
UF POQ7490 
UF PO7500 
UFPO75S10 
UFPO7520 
UFPO7530 
UF PO75S40 
UFPO7S50 
UFPO7S560 
UFPO7570 
UFPO7 S580 
UFPO7S30 
UF PO7600 
UFPO7610 
UF PO7620 
UFPO7630 
UFPO7640 
UFPO7650 
UF PO7660 
UFPO7670 
UF PO7680 
UFPO7690 
UFPO7700 
UEPO77 10 
UFPO7720 
UEPO7 730 
UFPO7740 
UEPO? 750 
UFPO7760 
UFPO7 770 
UFPO7780 
UFPO7790 
UF PO7800 
UFPO7810 
UFPO722z0O 
UFPO7830 
UFPO7840 
UFPO7350 
UFPO7860 
UFPO?78/70 
UFPO7880 
UFPO78390 
UFPO7900 
UFPO7910 
UFPO7920 
UFPO7930 
UFPO7940 
UFPO7950 
UFPO7960 
UFPO7970 
UFPO?72e0 
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300 
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100 


C «2 


COMMON / HARCOM / A(15,3),8(15,3),AMP(15,3),PH(15,3),NH 


DIMENSION THETA(80),VX(80) ,UA(60,9) ,RWK(9) 
N=1 
NR=1 


CALCULATE VELOCITIES AT FIVE BLAOE GRID RADII 


QO 300 M=3,9,2 4 
NR=NR+ 1 


CALCULATE VELOCITIES AT NSR ANGLES FOR EACH RADIUS 
00 200 KK=1,NSR 

FIELD POINT CENTERED ON BLADE GRIO 

CALL FPGRIO(X,Y.Z.N.M.KK, 12) 

IF (M.EQ.3) X1=XFP 

KK15=KK+15 

THETA(KK15)=TZFP 

RWK(NR)=RFP 

DTFP=360. /FLOAT(NSR) 


CALCULATE AXIAL VELOCITY FOR FIELO POINT AT RAOIUS RWK(NR) ANO 
ANGLE THETA(KK1S5). STORE THIS VELOCITY IN VX, 


CALL AMPCAL(VX(16),KK) 
NSR1IS=NSR+15 


STORE SMOOTHEO VELOCITIES IN NR-TH COLUMN OF UA 
CALL SMOOTH(NSR,NR,VX,THETA,UA) - 
CALCULATE NSR VELOCITIES AT PROPELLER HUB 
NRHUB = 1 

XFP=1.5*Xi 

RWK (NRHUB ) =RH- 

STEP THROUGH NSR ANGLES AT HUB 

DO 100 KK21,NSR 

KK 15=KK+15 
THETA(KK15 )=360. =FLOAT(KK-1)/FLOAT(NSR) 
TZFP=THETA(KK15) 

CALL AMPCAL(VX(16),KK) 

NSRIS=NSRt15 


NOW STORE SMOOTHED HUB VELOCITIES IN FIRST COLUMN OF UA 


acs ae 


UEPO7 330 
UF PO8000 
UFPO8010 
UFPO8020 
UF PO8030 
UFPO8O040 
UFPC8O50 
UFPO8060 
UFPO8070 
UF PO8O080 
UFPO8090 
UFPO8 100 
UFPO8110 
UFPO8 120 
UFPO8 130 
UFPO8 140 
UFPO8 150 
UFPO8 160 
UFPO8170 
UFPO8180 
UF PO8 1390 
UF PO8200 
UFPO8210 
UFPO8220 
UEPOS230 
UF PO8240 
UPFPOsS250 
UF PO8260 
UFPO8270 
UFPO8280 
UFPO8290 
UF PO8300 
UFPO8310 
UFPO8320 
UFPO8330 
UFPO8340 


- UFPO8350 


UFFPOBsGO 
UFPO8370 
UF PO8380 
UFPO8390 
UFPO8400 
UF PO8410 
UFPO8420 
UFPO8430 
UFPO8440 
UFPO8450 
UFPO8460 
UFPO8470 
UFPO8480 
UFPO8490 
UF PO85S00 
UFPO8510 
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500 


axe 


ann 


901 
302 
101 
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CALL SMOOTH(NSR,NRHUB, VX, THETA,UA) 


NOW CALCULATE VELOCITIES AT A NEXT RADII OUTSIOE THE PROPELLER 
DISK. "USE APP OF PROPELLER GRID CASE OF N21. M29. 


DIAGNOSTICS &8&666&&8&&88858h88 


GOTO 101 


CAVGUMPATE NSR VESCGTVIES AT t.1,1.4, AND 1.7 PROP RADII 
DG S00 KKK=1,7,3 

RFP=1.0+FLOAT(KKK)/10. 

NR=NR+ 1 

RWK(NR)=RFP 


CAECUEATE NSR VELOCITIES AT THESE RADII 


DO 400 KK=1,NSR 

KK15=KK+15 

THETA(KK15)=360. *FLOAT(KK-1)/FLOAT(NSR) 
TZFP=THETA(KK15) 

CALL AMPCAL(VX(16),KK) 


NOW STORE SMOOTHED VELOCITIES IN NR-TH COLUMN OF UA 
CALL SMOOTH(NSR,NR,VX, THETA,UA) 

WRITE(OUTPUT IN FORMAT FOR UNSTEADY PIEWAKE 
FORMAT(8F10.5) 

FORMAT({ 1615) 

WRITE(8,902) NR 

WRITE(8,901) (RWK(M),M=1,NR) 

WRITE(8,901) ((UA(KK,M),KK=1,NSR),M=1,NR) 

STOP FPUV UPON RETURN TO MAIN PROGRAM 

ISTOP=1 


RETURN Pa 
ENO 


UF PO8520 
UFR@aS530 
UFPO8S40 
UFPO8550 
UFPO8560 
UF PO8570 
UFPO8S8O 
UPFeOas 90 
UFPO8600 
UFPO8610 
UF PO8620 
UFPO8630 
UFPO8640 
UFPO8650 
UF PO8660 
UFPO8680 
UFPO8s IO 
UFPO8700 
UFPO8710 
UEPOS7 20 
UFPO8730 
UFPO8740 
UFPO8750 
UFPO8760 
UFPO8770 
UFPO8780 
UFPO8790 
UF PO8800 
UFPO8310 
UFPO&85&20 
UFPO8830 
UFPO8840 
UFPO8850 
UFPO&8860 
UFPO8870 
UFPO8880 
UFPO8890 
UFPO8S00 
UFPO8910 
UFPO8920 
UFPO8930 
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SUBROUTINE AMPCAL(VX,KK) 


weexeee eee anaeeeeaaeaeeczneaeaneeeaeeaanaeaeaeeeeegaanaeaaewkeeeneeeeaneseeaaneeneanuanne 


= x 
* SUBROUTINE AMPCAL WILL TAKE FIELO POINT GEOMETRY * 
* CALCULATED IN PIEOAT ANO CALCULATE THE INDUCEO * 
* TIME-AVERAGED AXIAL VELOCITY. IT WILL STORE THAT®* 
* VELOCITY IN VX(KK). * 
x Pd 
* NSR => NUMBER OF ANGLES PER REVOLUTION * 
* K => ANGLE INCREMENTS + 1 * 
* VX( KK) => WHERE IME-AVERAGEO AXIAL VEL STORED »* 
* AMP(1,1) => ZEROETH-HARMONIC AXIOAL VELOCITY = 
* COMPUTED BY FOUCAL * 
z = 
SERRATE RAERAER ABR EHTS HA REET SHENAE ARATE AEE SE TSS 


*** COMMON TO CONTAIN GEOMETRIC PARAMETERS 

COMMON / GEOCOM / MM,MN,NN,NW,NTORL,X(12,10),.Y(12,10),2(12,10), 
& XTIPC 10,11). YTIP( 10, 11), ZTIP( 10, 11),XH( 101.7). YH(101,7). 

& Z2H(101,7),XW(21,10),YW(21,10),Z2W(21,10),RH 


wae COMMON TO CONTAIN SINGULARITY STRENGTHS 


COMMON / SINCOM / GT(10,10,60),GB(10,9,60),GTW(10),S8(10,9),GTV 
& ,GAM(9,60),GMEAN(9),OKFW(20) ,HBWX,HBWY ,HBWZ 
s** COMMON TO CONTAIN FIELO POINT DATA, NBLAOE, IDENT, NSR, AND VSR 
COMMON / FPCOM / NBLAOE,NSR,XFP,RFP,TZFP,YBAR, ZBAR, 
& OBLAOE,OTFP,IDENT(18),TFP(60),VSR,IPLOT,OIAM 


«** COMMON TO CONTAIN VELOCITY OATA 
COMMON / VELCOM / VIX(60),VIT(60),.VIR(60).U(60,3),NSAMP 
x* 2 


COMMON TO CONTAIN HARMONICS COEFFICIENTS 


COMMON / HARCOM / A(15,3),8(15,3),AMP(15,3),PH(15,3),NH 


OIMENSION VX(1) 
sae RUN THROUGH ANO FPSTEP FOR EACH TIME STEP 
00 100 K=1,NSR 
wee FPSTEP CALCULATES VELOCITY INOUCED AT A FIELO POINT BY ONE 
wae BLAOE AT TIME STEP K 
100 CALL FPSTEP(K) 


5 BS 


UFPO8940 
UFPOSSSC 
UFPO8S60 
UFPO8970 
UPPOs sso 
UFPOS9SSO 
UFPOSO00 
UFPOSOIG 
JEPCIC 20 
UFPOSO30 
UFPOSO40 
UFPOSOSO 
UFPOSO6O 
UFPOSO70 
HEPOSOS©e 
UFPOSOSO 
UE POS 100 
UFPOS110 
URPOS1Z¢ 
UF POS 130 
UFPO9S140 
UPPED 1S0 
UFPOS 160 
UFPOS170 
UFPOS 180 
VErOsS1S9 
UFPOS200 
UFPOSZ1O 
UFPOS220 
UFPOS230 
UFPO9S240 
UE PCSZ5© 
UERPOSZ5© 
UFPOS270 
UFPOS280 
UF POS2390 
UF POS200 
UF POSS1O 
UFPOS320 
UFPO9S330 
UFPOS340 
UFPOS3S50 
UFPOS360 
UFPOS9370 
UF POS Seo 
UFPOS330 
UFPOS400 
UFPO9410 
EPCS 470 
UF POS430 
UFPOS440 
UFPOS450 
UFPOS460 


. 
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SUMVEL SUMS FIELD POINT VELOCITIES INOUCED BY EACH BLADE 


CALL SUMVEL 


FOUCAL GENERATES HARMONIC COEFFICIENTS FOR INOUCED FIELD POINT 


VELOCITY AS PROPELLER ROTATES. 


CALL FOUCAL 


STORE TIME-AVERAGED (ZERO-HARMONIC) VELOCITY IN VX 


VX(KK)=AMP(1,1) 
RETURN 
END 
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UFPO9470 
UFPOS480 
UFPOS490 
UF POSS00 
UFPO3S510 
UFPO9S20 
UFFOSs20 
UFPOS540 
UFPOSSS5O0 
UFPOS560 
UFPOSS70 
UFPOSS80 
UEPOSS3O 
UF POSE00 
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100 
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200 


C «2 


300 
201 
=)OPs 


C «a 


SUBROUTINE SMOOTH(NSR.NR, VX, THETA,UA) 


oe oe oe oe eo ee oe oe ee oe oe oe ee oe oe oe oe oe ee ee ee ee eo ee oe ee ee ee ee ee ee oe oe oe me 


x 
= 
= 
* 
" 
® 
x 
ma 
x 
® 
® 
* 
4 
x 
® 
* 
= 
® 
= 
® 
x 
x 


DIMENSION VX(80),UA(60,9),THETA(80),TH(60),AVX(250) 
COPY LAST 15 STEPS IN REV INTO FIRST 15 ENTRIES IN VX 


00 100 K2=1,15 


x 


SUBROUTINE SMOOTH REAOS INOUCED VELOCITIES AT ONE® 


RAOTUS FROM PIEDAT. 


IT CALCULATES BY 


INTERPOLATION A SERIES OF NSR EVENLY SPACEO 


AXIAL INOUCED VELOCITIES SUITABLE FOR INPUT INTO 


UNSTEADY PIEWAKE. 


NSR 
NR 


VXx(1 - 5S) 


Vx(6 - NSR+5) 
VX(NSR+6 


THETA(K) 
UA(K,M) 
TH(K) 
NIN 
KPOS 


=> 
=> 


*=> RAW VELOCITIES 5 STEPS BEFORE ZERO 


=> 
=> 
=> 
=> 
=> 


KSTART,KSTOP 


AVX 


KNSR=K+NSR 
THETA(K)=THETA(KNSR)-360. 


VX(K)=VX(KNSR) 
COPY FIRST S TIME STEPS IN REV INTO S ENTRIES AFTER REVOLUTION 


00 200 K2=16,20 


KNSR=K+NSR 
THETA(KNSR)=THETA(K) + 360. 


VX(KNSR)=VX(K) 


=> 


NUMBER OF ANGLES PER REVOLUTION 
RAOIUS NUMBER 1=HUB, ETC. 


BLADE ANGLE. 


* 


¢ #&# * # # © © & 


=> RAW VELOCITIES IN ONE REVOLUTION® 


- NSR+20) => RAW VELOCITIES 15 TIME 


TIME STEPS AFTER ONE REVOLUTION 
ANGLE OF VELOCITY vx(K) 

SMOOTH VEL AT RAOIUS M, ANGLE K 
EVENLY SPACED ANGLES OF UA(K,M) 
NUMBER OF ANGLES FED TO UGLYOK 

K OF FIRST POS. ANGLE IN THETA(K) 


=> FIRST ANO LAST INOEX OF THETA(K) 


AND VX(K) SENT TO UGLYOK 
SPUINES CURVE COcFFICIENTS OF VX 


R2eeeeaecgceaneeaanecneeeeeaneeeacnewe eee acaaeeeeceeneneaneesteswanwe ewe ene 


SET UP EVENLY SPACED ANGLES 


DO 300 K=1.NSR 


TH(K)=360. =FLOAT(K-1)/FLOAT(NSR) 


FORMAT(8F 10.5) 
FORMAT(161S) 


FIND FIRST POSITIVE ANGLE, SAVE INDEX 


ee 


* © # © # & © &# © © © & 


ANO THETA 


UFPO96 10 
UF PO9620 
UFPOS630 
UFPO9S640 
UF POS650 
GProy ogo 
UF PO9670 
UFPO9S680 
URE OVE90 
UFPO9700 
UFPO9710 
UFPO9720 
UFPOS730 
UFPO9740 
UF PO9750 
UEPOS 760 
URPO3770 
UFPOS780 
CEPGS790 
UF POS 800° 
UFPOS810 
URPOSs 20 
UFFPCSE3¢O 
UFPOS840 
UF PO9850 
UFPO9S860 
UFPO9870 
UF PO9880 
URPOs230 
UFPOS300 
UFPO99 10 
UFPOSS20 
UF PO9930 
UFPO9S940 
UFFPOSSs0 
UFPOS960 
UFPOS3970 
UFFOSSSO 
UFPOSS30 
UFP 10000 
UFP1O010 
UFP 10020 
UFP 10030 
UFP 10040 
UFP 10050 
UFP 10060 
UFP 10070 
UFP 10080 
UFP 10090 
UFP10100 
UEP 1OT1©0 
UFPTO120 
UFP 10130 





Cc 


Cc 


400 
300 


zu 


00 400 K=1,NSR 
KPOS=K 

IF (THETA(K).GT.0O.) GOTO SOO 
CONTINUE 

KSTART=KPOS- 1 

NSR2=NSR/2 

KSTOP=KPOS+NSR2+1 
NIN=KSTOP-KSTART+1 


FINE SPLINE COEFFICIENTS ANDO SMOOTH VEL FOR FIRST 180 DEGREES 


CALL UGLYOK(NIN, 1,1, THETAC(KSTART) ,VX(KSTART),O.,0.,AVX) 
CALL EVALOK(NIN,NSR2, THETAC(KSTART),TH(1),UA(1,NR),AVX) 


FIND SPLINE COEFFICIENTS AND SMOOTH VEL FOR LAST 180 OEGREES 


KSTART=KSTOP-2 

KSTOP=KSTART+NIN- 1 

NSR2P=NSR2+1 

CALL UGLYOK(NIN, 1,1, THETAC(KSTART) ,VX(KSTART),O.,0. ,AVX) 

CALL EVALDK(NIN,NSR2, THETA(KSTART),TH(NSR2P) ,UA(NSR2P,NR),AVX) 
RETURN 

END 
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UEP 16 140 
UFP 10150 
UFP 10160 
UFP 10170 
UFP 10180 
UFP 10190 
UFP 10200 
UFP 10210 
URPIO226 
UEP T1C23¢C 
UFP 10240 
UFP 10250 
UFP 10260 
UFP 10270 
UFP 10280 
URE 10290 
UFP 10300 
UFP1O310 
NFP 1C3s26 
UFP 10330 
UFP 10340 
UEP TOS oo 
UFP 10360 
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. PIEWAKE Program Listing 
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110 
103 
aS 
114 
13 
111 
210 
200 


C «= 


THIS PROGRAM COMPUTES EFFECTIVE WAKES BY ONV METHOO 


Save eepeeeseeee keeezazrwavaeswewseeeaeeeeaveaeneeeweaea caer weae ec eewaoeweeavrvrananee 


= 


* *£ © # @ © £ © © © © &# #©# © © © © © # #8 8 8H © # 8h He hUthlhUhthlUrthlUhhlUhlhlUrtllU 


THE ORIGINAL PIEWAKE WAS WRITTEN BY R. 
1981. 
VERITAS METHOO OF COMPUTING EFFECTIVE WAKE. 
PROGRAM HAS BEEN MOOIFIEC BY R. 


IN ABOUT 


= 
VAN HOUTEN# 
IT 1S BASEO ON THE OET NORSKE . 
(Enpes 
JAMISON IN APRIL,* 


1982, TO ACCOMOOATE UNSTEAOY PROPELLER-INDUCEO * 


AXIAL VELOCITIES FROM UFPV. BASICALLY, BOTH * 
PROGRAMS OIVIOE THE WAKE INTO A NUMBER OF PIE- * 
SHAPEO SEGMENTS. EACH SEGMENT IS ASSUMEO TO = 
CONTRACT IN AN AXISYMMETRIC MANNER AS GIVEN BY * 
T. HUANG. * 
NP => NUMBER OF PIE-SHAPEO SEGMENTS * 
NX => NUMBER OF NOMINAL WAKE RAOII * 
RX(M) => NOMINAL WAKE RAOIUS M (UP TO 11) * 
VO(K,M) => NOMINAL AXIAL VELOCITY, SEG K, RAO M* 
VT(K,M) => NOM TANG VEL, SEG K, RAO M * 
VR(K,M) => NOM RAO VEL, SEG K, RAO M * 
VTE,VRE => VT & VR INTERPOLATEO AT EFF WAKE RAQ* 
UASEG(M) => UA IN A PIE-SEGMENT. ALSO USEOQ TO -* 

STORE A PIE -SEGMENT OF VT s 
UX(M) => VO IN A PIE-SEGMENT. ALSO USEOQ TO * 

STORE A PIE-SEGMENT OF VR * 
AUX, AUASEG=> SPLINE COEFF. FOR VR,VT, RESP. * 
RE(M) => RAOII OF EFFECTIVE WAKE VELOCITIES * 
THETA(K) => ANGLE IN OEGREES OF PIE SEG K * 
NA => NUMBER OF PROP-INO, AXIAL VEL RAOII * 
RA(M) => PROP-INO AXIAL VEL RAOIUS M * 
UA(K,M) => PROP INO AXIAL VEL, SEG K, RAO M = 
VE(K,M) => EFFECTIVE WAKE VEL, SEG K, RAD M = 
VL => AVERAGE VELOCITY * 
VLAVN => VOLUMETRIC AVERAGE NOMINAL VELOCITY * 
VLAVE => VOLUMETRIC AVERAGE EFFECTIVE RSet TY 
BLOCK => EFFECTIVE BLOCKAGE 


2 


aeeseeaeezeaee eee ee eeaeeeeeezeeecnagngcaneaeezweaee eee eseeeweneseaeeeeeeeeaeaeesee 


OIMENSION VO(60,11),RX(11),UX(60),RA(11),.UA(60, 11). UASEG(60) 
VR(60,11),VT(60, 11), AUASEG(44),AUX(44),VTE(11) 
DIMENSION THETA(60).,VE(60.11)., IOENT(18),UERX(60), VRE( 11) 


DIMENSION RE(11), 


COMMON RX,UX.RA,UASEG, UERX,RE 


FORMAT(S8F 10.5) 
FORMAT (/’ 
FORMAT(/’ 
FORMAT(/’ 
FORMAT (/’ 


FORMAT(18A4) 
FORMAT( 1615S) 


FORMAT (8F 10.5) 


REAQ IN NOMINAL WAKE OATA. 


AVERAGE VELOCITY *’,F5.3) 
VOLUMETRIC AVERAGE NOMINAL VELOCITY =’,F5.3) 
VOLUMETRIC AVERAGE EFFECTIVE VELOCITY =’,F5.3) 
EFFECTIVE BLOCKAGE *’,F7.3, 


PERCENT’ ) 


= 6 = 


SAVE RAQIAL ANO TANGENTIAL 


PIT EOCC1S 
PITEQCOZS 
PTEOOO30 
PTEQOO4O 
PIECOGSOG 
PIEQOO6O 
PI EQOO70 
PTEQOO80 
PTECCOIO 
PTEOO100 
PIEO0O110 
PIECE IZE 
PIEOO 130 
PIEO00140 
PTEQOO150 
PIEGOT6O 
PITECOT7 © 
PTEOO180 
PILEQO1S0 
PTEOQO200 
PIEOOZ IC 
PITEQO0220 
PIEGO2Zs0 
PIEQOO240 
PIEGO2ZS0 
PIEQO260 
PIEOOZ70 
PTEOO280 
PLECOZSO 
PT EOO300 
PIEQ0310 
PIEOQO320 
PIEQOO330 
PITEOQO340 
PIEQO33SC 
PTEGOS6® 
PIECO27C 
PTEOO380 
PITEOO2S0 
PI EQO400 
PTEOO410 
PIEQ0420 
PIEOO430 
PIEOO440 
PIEQOO4SO 
PIEOO460 
PIEQO470 
PIEOO480 
PIT EQOO490 
PTEOOSOO 
PIECO2 10 
PL ECOGs20 
PT EOQOS30 





c VELOCITIES. THEY ARE NOT NEEDEO FOR EFFECTIVE WAKE PIEOQOS40 
Gc CALCULATIONS BUT WILL BE ECHOED FOR WKPRDC. PIEOOSSO 
c PIEOOS6O 
REAO(10,111) (IOENT(I),I=1,18) PIEOO570 

REAO( 10,210) NX PIEOOS8O 

REAO( 10,210) NP PIEOOS390 

READ( 10,200) (RX(M) ,M=1,NX) PIEQO600 

DO 20 M=1,NX PIEOQO610 
REAO(10,200)((THETA(K),VO(K,M),VT(K,M),.VR(K,M)),K=1,NP) PIEQO620 

20 CONTINUE a PIEOCO630 

S é PIEOQO640 
C *=#* REAO IN PROPELLER-INOUCED VELOCITIES FROM EITHER FPUV PIEOO6S5O 
C OR A SEPARATE OATA FILE PIEQO6E6O 
@ PIEOO670 
900 FORMAT(/’ 00 YOU WISH TOD USE STEAOY OR UNSTEAOY VELOCITY PROFILE? ’PIEOO680 
& /’ ENTER 1 FOR STEADY, 2 FOR UNSTEADY’ ) PIEQO690 

8 WRITE(6,900) PIEQO700 
READ(S,*) ISTEO PIEQO710 

IF (ISTED.GT.2.O0R.ISTEO.LT.1) GOTD 8 PIEOQO720 
REAO(8,210) NA PIEOO730 
REAO(S.200) (RA(I).12=1.NA) PIEOO740 
RH=RA(1) PIEOQO7S5SO 
RPROP=1. PIEOO760 
lFursneOleG@.1) GOTD 22 PIEOO770 
REAO(8,200) ((UA(KK,M),KK=1,NP).M=1,NA) . PIEOO780 

GDTO 92 PIEOO790 

22 REAO(8,200) (UA(1,M),M=1,NA) PIEOO800 
DO 21 KK=2,NP PIEOO810 

“ OO 21 M=1,NA : PIEOO820 
21 UA(CKK,M)=UA(1,M) PIEOCO830 

c PIEOQO840 
C #*2 EXTRAPOLATE NOMINAL WAKE DATA TO: HUB PIEOO8SO 
@ PIEOO860 
92 WRITE(6, 100) PIEOO870 
100 FORMAT(/’ OO YOU WISH TO EXTRAPOLATE VELOCITY DATA TD THE HUB?’°/ PIEOO880 

1 %’ ENTER 2 FOR NO, 1 FOR LINEAR, O FDR CONSTANT’ ) PIEOCO830 
REAO(S,*)IXTRAP PIEOQOS90O 

NX 1=NX PIEOCO910 

IFC IXTRAP.GT.1)GO TO 28 PIEOQO920 

c PIEOO930 
C **#* LINEAR EXTRAPOLATION OF NOMINAL WAKE OATA TO HUB. PIEQO940 
Cc ; PIEOO950 
NX=NX+ 14 PIEOO960 

OD 25 K=1,NP PIEOOS70 

OO 23 M=1,NX1 PZEOO980 
MR=NX-M PIEOO990 

MRP 1=MR+ 14 PIEOQ1000 
VT(K,MRP1)=VTC(K.MR) PIEO1010 
VR(K,;MRP1)=VR(K,MR) PIEO1020 

23 VO(K,MRP1)=VO(K,MR) PIEO1030 

25 VO(K,1)=VO(K,2)+IXTRAP*(VO(K,2)-VO(K,3))/(RX(1)-RX(2))*(RH-RX(1)) PIEO1040 

DO 26 M=1,NX1 PIEO10S0O 
MR=NX-M PIEO1060 


= ons 





a00 


a00 


aan 


o0a0 


anngn 


a00 


26 


28 


2a 


30 


Ce 


60 


70 


222 


40 
30 


MRP 1=MR+ 1 

RX (MRP 1)=RX(MR) 
RX(1)=RH 

VL=O0. 

VLAVN=O. 
VLAVE-O. 


CYCLE THROUGH EACH PIE SEGMENT 


DO 50 K=1,NP 
DO 30 M=1,NX 
UASEG(M)=UA(K,M) 
UX(M)=VO(K,M) 


WKMOOD PERFORMS THE HUANG MODIFICATION ON ONE PIE-SEGMENT 


CALL WKMOD(NX,NX,NA,1.RPROP,VOL,VOLAVN, VOLAVE ) 
INTERPOLATE VT AND VR AT EFFECTIVE WAKE RADII 


DO 60 M=1,NX 
UASEG(M)=VT(K,M) 

UX(M)=VR(K,M) 

CALL UGLYDK(NX,1,1,RX,UASEG,O..0O. , AUASEG) 
GALL UGUYDRONX. 1. 1.RX.UX.0..0.,AUX) 

CALL EVALDK(NX,NX,RX,RE,VTE,AUASEG) 

CALL EVALDK(NX,NX.RX, RE, VRE, AUX) 

DO 70 M=1,NX 

VT(K M)=VTE(M) 

VR(K,M)=VRE(M) 


ACCUMULATE VOLUMETRIC AVERAGE DATA 


VL=VL+VOL : 
VLAVN=VLAVN+VOLAVN 

VLAVE =VLAVE+VOLAVE 

DO 40 M=1,NX 

VE(K,M) =UERX(M) : 

CONTINUE S 


COMPUTE VOLUMETRIC AVERAGES 


VL=VL/NP 

VLAVN=VLAVN/NP 
VLAVE=VLAVE/NP 
BLOCK=(1.-VLAVE/VLAVN) #100. 


WRITE EFFECTIVE WAKE DATA IN FILE 11 TO BE PROCESSED BY WKPROC 
WRide@ii. 111) (IDENT(I),2=1, 18) 
WRITE( 11,210) NX 


WRITE( 11,210) (NP,M=1,NX) 
WRITE(11,200)(RE(M),M=1,NX) 


aoe) 


PIEO1O70 
PLEO ICSC 
PIEO1030 
PITEO1100 
PIEOT110 


“PIEO1120 


PIEO1130 
PIEQ1140 
PIEQ1150 
PIEQ1160 
PIEQ1170 
PIEO1180 
PIEO1190 
PIEO1200 
PIEOQ1210 
PIEO1220 
PIEQ1230 
PIEO1240 
PIEQ1250 
PIEO1260 
PIEO1270 
PIEO1280 
PIEO1290 
PIEO1300 
PIEO1310 
PIE01320 
PIEOQ1330 
PIEO1340 
PIEO1350 
PIEO1360 
PIEO1370 
PIEQ1380 
PIEO1390 
PIEO1400 
PIEO1410 
PIEQ1420 
PIEOQ1430 
PIEOQ1440 
PIEO1450 
PIEO1460 
PIEO1470 
PIEQ1480 
PIEO1490 
P1EO1500 
PIrEGis1¢ 
PIEO1520 
PIEQ1530 
PIEO1540 
PIEQ1550 
PIEO1S60 
PIEQ1570 
PIEO1580 
PIEQ1590 





aan”n 


10 


DO 120 M=1,NX 


Whe @nl!, 2co)((THETA(K), VEC(K M) (VICK, M)VR(K,M)),K=1,.NP) 


CONTINUE 
WRITE AVERAGES TO TERMINAL 


WRITE(6,103)VL 
WRITE(6,.113)VLAVN 
WRITE(6,114)VLAVE 
WRITE(6,115)8LOCK 
STOP 

END 


+ GE) 


PIEO1600 
PIEQIS1¢ 
PIEO1620 
PIEQ16E30 
PIEO1640 
PIEGRS Ie 
PIEOQ1660 
PIEOQ1670 
PIEO16Eao 
PIEO?6E90 
PIEOQ1700 
PIEO47 10 
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SUBROUTINE WKMOO(N,NX,NA,IOW,RPROP, VOL, VOLAVN, VOLAVE) 


ZQeeeaeeaeveeneaeeegeavwveve ve eaeveeceervee ee eevee RB aeweeweveazeaeweeneczcenenseveevs & 


SUBROUTINE WKMOO CALCULATES THE THOMAS HUANG 
CONTRACTION OF THE NOMINAL ANO INOUCEO WAKES 
PASSEO TO IT THROUGH COMMON. IT ALSO COMPUTES 
VOLUMETRIC AVERAGES. 


NX => NUMBER OF NOMINAL RADOTI 

RH => RADIUS OF HUB 

RT => RADIUS OF OUTERMOST NOMINAL VELOCITY 
UX, UXR => NOMINAL WAKE VELOCITIES 

UE, UEAR => EFFECTIVE WAKE VELOCITIES 

UP => APPARENT WAKE = EFF. WAKE + INO. WAK 
UA,UAR => INOUCEOD VELOCITIES 

IOw = 1 => NO TUNNEL CORRECTIONS 

Lier => OPEN WATER ITERATION INOEX 

ITERAT => HUANG CORRECTION FACTOR INOEX 

SG, 0. F => HUANG CORRECTION FACTORS. SEE PAPER 


AUA,AUX,AUE=> SPLINE COEFF. FOR INOUCEO, NOMINAL, 
AND EFFECTIVE VELOCITY FIELOS. RESP 


* *# © © @# © © © © © © © #@ © # © © © 8 * 8 oe hele 
¢* ¢ © # # # # # # © © MM © © @ © @ © © © @& © & 


RPROP => RADIUS OF PROPELLER 
RE(M) => SPECIFIED EFFECTIVE WAKE RADII 

VOL => AVERAGE VELOCITY 

VLAVN => VOLUMETRIC AVERAGE NOMINAL VELOCITY 
VLAVE => VOLUMETRIC AVERAGE EFFECTIVE VELOCITY 


* 


Raeswsenrneveeeveeenengveeeveageaezce eaeveevenanesewveevevaeeecanekawereoeezk een eseveveveeee 


OIMENSION RX(11),R(011),RA( 11), UX(60),UXR(60),UA(60) ,UP(60),UE(60) 


OIMENSION RE(11),UAR(60),RP( 11) 
DIMENSION AUX(240),AUA(240) 
OIMENSION AUE(240),UERX(60) 
COMMON RX,UX,RA,UA,UERX,RE 


USE R(I) FOR NOMINAL RAOII 


N=NX 
Ni=N-1 
RH=RX(1) 
RT=RX(NX) 


a 


INTERPOLATE TO FINO NOMINAL VELOCITIES AT NOMINAL RADII 
CALL UGLYOK(NX,1,1,RX,UX,0O.,0O. , AUX) 

DO 10 I=1,N 

R(1)=(RT-RH)/(N-1)*#(1-1)+RH 

R(I)=RX(1) 

CALL EVALOK(NX,N,RX,R,UXR,AUX) 

INTERPOLATE TO FINO INOUCEOD VELOCITIES AT NOMINAL RAOII 


CALI UGEYDK(NA,1,.1.RA.UA.O..0. ,AUA) 


> 9h - 
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PLEQ17 30 
PTEO1740 
PLEG1750 
PTEO1760 
PIEOI770 
PIEOQ1780 
PLEGiI7s°O 
PTEGTSCO 
PTEO18 10 
PIEG1S20 
PIEO1ZS0 
PIEOQ1840 
PIEOle2¢ 
PIEOQIS260 
PTEOIBTO 
PIEOQI8E0O 
PIEO1890 
PIEO19Ce 
PIEO1910 
PIEOQ1920 
PTEOIS30 
PIEO1940 
PLEO lZo° 
PIECISSO 
PIEOQ1970 
PIE@G1 S80 
PIEO1990 
PIEOZ000 
PIEO2010 
PIeO2Z2020 
PIEG2Z030 
PTEO2040 
PIEO2050 
PITEO2060 
PITEQO2070 
PTEOCZO0E0O 
PIECZOSe 
PIEO@Z100 
PITEOQ2110 
PIEOZ120 
PTEOZ1320 
PTEOQ2140 
PIEOZ15¢C 
PIEOZ16C 
PIEOZ1I7TC 
PLEOZ1SO 
PIEOGZ120 
PIEO2200 
PIEGZ 210 
PILEGZ220 
PABOZ230 
PIEO2240 
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SO 


CALL EVALOK(NA,N,RA,R,UAR,AUA) 


FORMAT(8F 10.5) 


OUTSIOE PROPELLER OISK, NOMINAL WAKE 


IUPDAT=0 

ITER=1 

ITERAT=1 
UE(N)=UXR(N) 
UP(N)=VE(N)+UAR(N) 


ESTIMATE EFF. WAKE MARCHING FROM OUTSIOE TO HUB 


DO 30 IT=1,N1 
K=N-I 


ERRECTIVE WAKE 


PIEO2230 
PIEO2260 
PIEQ2270 
PIEO2280 
PIEO2290 
PIEO2300 
PIL EG2Z310 
PIEOZ320 
PIREOZs20 
PIEQ2340 
PITEO2Z350 
PIEO2360 
PIECZ37©¢ 
PILEG2Z23©0 
PIEO2390 


UE(K)=SORT( (UE(K+1)+(UAR(K+1)+UAR(K))/2. )=*2+UXR(K) =#*2-UXR(K+1)==*2PI1E02400 


1)-(UAR(K+1)+UAR(K))/2. 
UP(K)=UE(K)+UAR(K) 
CONTINUE 


COMPUTE HUANG’S CORRECTION FACTORS 


RP(1)=RH 

VOL =O. 

00 40 K=1,N1 
B=2*UP(K+1)+UP(K) 
C=-RP(K)=(UP(K+1)-UP(K) ) 


PIEQO2410 
PIEO24:2@ 
PIEO2430 
PIEOQ2440 
PIEOQ2450 
PIEQ2460 
PIEOQ2470 
PIEO2480 
PIEO2490 
PIEG2300 
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F=(R(K+1)**2-R(K)*#2)*(2"UXR(K+1)+UXR(K))-R(K)*(UXR(K+1) -UXR(K)) #*(PIEO2520 


1R(K+1)-R(K)) 
VOL =VOL+F 


O=-RP(K)*=*22(UP(K+1)+2*UP(K))-F 


FINO RAOIUS FOR (K+1)ST APPARENT WAKE 


IF ((C*C-4*8*D).LT.O.) GOTO 80 


RP(K+1)=(-C+SORT(C#*2-4*B8+0))/2./B 
. IFC ABS(RP(N)-RT).LT.(.0001*RT))GO TO 60 


IUPOAT=0 


TOW=O0 => TUNNEL CORRECTIONS 


x1 => OPEN-WATER (I.E., 


1F(l@W.EBO.1) GO TO 60 
IF (ITER.NE.1)GO TO SO 
ITER=2 

UENOLO=VE(N) 
RPNOLD=RP(N) 
UE(N)=UE(N)*®.98 

GO TO 20 

CONTINUE 

ITER=ITER+1 
IF(ITER.GT.10)GO TO 80 
UEN=UE(N) 


INFINITE FLUID) 


= Oa 
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PITEO2740 
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UE(N)=UE(N)+(RT-RP(N) )/(RP(N)-RPNOLD) *(UE(N)-UENDLD) 
RPNDLO=RP(N) 

UENDLD=UEN 

GD Gan 20 


ITERATE ON HUANG CORRECTION FACTORS, 3 TIMES FOR OPEN WATER 
IF(IUPOAT.EQ.1)GO TO 70 
IF(IOW.EQ.1.AND.ITERAT.EQ.5)GD TO 70 

ITER=1 

ITERAT=ITERAT+1 

IFC ITERAT.GT.10)GO TO 890 

IUPOAT=1 

FIND INOUCEOD VELOCITIES AT APPARENT WAKE RADII 
CALL EVALOK(NA,N,RA,RP,UAR,AUA) 

GD jO 20 

CONT INUE 

CALCULATE EFFECTIVE WAKE RADII 

NXM1=NX- 4 

DO 100 M=1,NX 

RE(M)=RH+(RPRDP-RH) *FLOAT(M-1)/FLDAT(NXM1) 


IF OUTERMOST APPARENT RADIUS IS INSIDE PROPELLER DISK, THEN 
LINEARLY EXTRAPOLATE TD PROPELLER RADIUS 


IF (RP(N).GE.RPROP) GOTO 110 

UE(N)=UE(N1) + (UE(N)-UE(N1)) *(RPROP-R(N1))/(RP(N)-RP(N1) ) 
RP(N)=RPROP 

INTERPOLATE EFFECTIVE WAKE VELOCITIES AT SPECIFIED RADII 


GAUI UGEYDK(N., 1.47.RP.UE.0..0..AUE) 
CALL EVALOK(N.NX.RP,RE,UERX, AUE ) 


COMPUTE VOL = AVERAGE VELOCITY 
VDOL=VOL/3./(RT**2-RH=*2) 
COMPUTE VOLUMETRIC AVERAGE NOMINAL VELOCITY 


CALL INTEOK(NX,RX,RH,RPROP,YOX.XYOX,XXYOX, AUX) 
VDLAVN=XYOX *#2/(RPROP * #2-RH*22) 


COMPUTE VDLUMETRIC AVERAGE EFFECTIVE VELOCITY 
CALL INTEOK(N,RP,RH,RPROP,YOX,XYDOX,XXYDX, AUE ) 


VOLAVE=XYOX*2./(RPROP=*2-RH=2) 
RETURN 
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C 
C 


80 


TUNNEL CORRECTION OR HUANG CORRECTION ITERATIONS FAIL TO CONVERGE PIEO3310 


STOP 
END 
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